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Preface 


Theses  in  the  past  developed  flight  controllers  for 
rigid  body  aircraft  .  In  this  thesis  one  digital  flight 
controller  is  developed  for  both  a  rigid  body  and  flexible- 
fighter  aircraft.  That  is  a  bending  mode  is  combined  with 
the  dynamic  of  the  rigid  body  aircraft  to  form  the 
flexible-fighter . 

I  would  like  to  thank  both  Dr  Robert  Schwanz  and 
Mr  Finley  Barfield  of  the  Flight  Dynamics  Laboratory  for 
their  assistance  and  guidance  in  the  development  of  the 
bending  mode  equation  for  the  flexible-fighter  aircraft. 
Without  their  help  and  suggestions,  the  development  of  the 
flexible-fighter  aircraft  would  not  have  been  possible. 

The  technique  used  in  the  design  of  the  flight 
controller  is  based  on  the  theory  developed  by  Professor 
Brian  Porter  from  the  University  of  Salford,  England.  I 
would  like  to  express  my  appreciation  to  my  thesis  advisor 
Dr  John  J.  D'Azzo  and  to  my  reader,  Dr  Constantine  H.  Houpis, 
both  of  whom  are  Professors  of  the  Department  of  Electrical 
Engineering  at  the  Air  Force  Institute  of  Technology. 

This  thesis  was  accomplished  in  parallel  with  four 
fellow  students.  I  would  especially  like  to  thank  Capt. 
William  Locken  who  also  helped  in  the  development  of  the 
flexible-fighter  aircraft  model.  Also,  my  sincere  thanks 
to  Lt.  Jeff  Simmers,  Lt.  Brian  Mayhew,  and  Lt.  Scott  Feldmann 
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who,  along  with  Capt.  Locken,  comprised  the  discussion 
group  which  aided  in  accomplishing  my  objectives. 

Finally,  but  most  of  all,  I  want  to  express  my 
love  and  appreciation  to  my  parents  and  the  rest  of 
the  family  who  supported  me  during  my  stay  at  AFIT. 

Lt  Marc  L.  Hoffman 
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8  Pitch  angle 

8C  Pitch  angle  command 

0  Pitch  angle  error 

T(X)  Asymptotic  transfer  function  matrix 
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Abstract 

-p  This  thesis  evaluates  the  effect  of  a  digital  flight 
controller  on  a  flexible-fighter  aircraft  (FFA)  at  <ihe 
flight  condi tion-of^O .8  Mach  at  30,000  feel.  The  two 
maneuvers  incorporated  in  this  flight  controller  are  the 
g-command  pull-up  maneuver  and  the  pitch-pointing 
maneuver^.  Using  a  control  law  that  is  based  on  the  control 
system  theory  developed  by  ^refeseor^Brian  Porter  ^t  the 

^  '  "N 

University  of  Salford,  England, >a  digital  flight  controller 

I  y‘:  J 

is  first  designed  for  a  rigid  body  aircraft,,  The -same 

„ —  -  -  -  -  _-■■■■■  -  i thsrs 

design  ^.s  then  used  on  a*!  f lexible^f ightei  ~alTcraTlr^/^Cfte^ - 

v  ^  A,  , 

—  -f,  effect  of  the  controller  for  the  rigid  body  design  when 
used  on  the  f  lexible-fighter  <are~avaluated 
^  >  The  results  showed  excellent  performance  by  the  flexible-  > ^  ^  ^  • 

r  - 

-f±ghtSr  aircraft 

The  flexible-fighter  model  is  based  on  the  rigid  body 
dynamics  of  the  AFTI/F-16  combined  with  the  first  bending 
mode  of  a  flexible  aircraft.  The  bending  mode  of  the  F  L  k 
-flexible  -aircraft  was  selected  to  correspond  to  expected 
characteristics  in  a  typical  fighter  aircraft. 

Incorporating  a  time  delay  in  the  controller  to 
represent  the  expected  microprocessor  computational  time 
delay  did  not  significantly  affect  the  simulation  time 


xx 


F  r  A  •  ' 

response  of  either  the  rigid  body  or  the  fiexifete-fighter 
aircraft.  '  +  x 

To  ^determine  the  acceptance  of  the  design  by  a  pilot* 
the  Neal-Smith  Criterion  is  used.  The  results  of  ^ 

applying  this  criterion  showed  that  the  pilot  would  accept 
the  design  for  both  the  rigid  body  and  the  flexible-fighter 
aircraft. 
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EVALUATION  OF  A  DIGITAL  FLIGHT  CONTROLLER 
FOR  A  FLEXIBLE-FIGHTER  AIRCRAFT 


I .  Introduction 


Background 


Uodern  flight  controller  designs  involve  strongly 
coupled  Multiple  Input-Multiple  Output  (MIMO)  non  linear 
plants  with  large  parameter  variations  due  to  a  wide 
range  of  Mach  numbers  and  altitudes.  Many  controllers  in 
the  past  have  been  designed  on  the  basis  of  Single  Input- 
Single  Output  (SISO)  system  design  using  conventional 
techniques  such  as  root-locus  and  frequency  response 
analysis  (Ref.  1).  The  state  and  output  equations  for  a 
SISO  system  have  the  form 


x  *  Ax  +  bu 


where 


y  -  Cx 

• 

x  is  an  n  x  1  column  vector, 

A  is  an  n  x  n  plant  coefficient  matrix, 

x  is  an  n  x  1  state  vector, 

b  is  an  n  x  1  control  matrix, 

u  is  a  one  dimensional  control  vector, 
£  is  a  one  dimensional  output  vector, 

C  is  an  1  x  n  output  matrix, 


(1.1) 


(1.2) 
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n  is  the  set  of  first-order  differential 
equations  with  n  independent  states. 

The  SISO  system  design  methods  have  been  around  for  sometime 
and  many  examples  and  short  cuts  have  been  developed  for  the 
different  design  techniques.  It  can,  therefore,  be  used  by 
a  control  designer  by  replacing  the  MIMO  system  with  a 
number  of  separate  SISO  systems  (Ref.  2).  Because  the  SISO 
systems  are  intercoupled ,  this  approach  requires  an  iterative 
approach  which  can  make  the  design  process  tedious. 

Because  of  the  advancement  in  computational  technology, 
design  procedures  for  MIMO  systems  can  be  more  readily  applied. 
For  this  reason,  staying  with  a  MIMO  system  design  approach 
offers  a  greater  amount  of  control  and  greater  precision  for 
the  design  of  Digital  Flight  Control  Systems  (DFCS).  Many 
different  methods  and  theories  have  been  developed  for  multi- 
variable  control  designs.  One  recent  method  has  been 
developed  by  a  Professor  Porter  at  the  University  of  Salford, 
England,  exploiting  high-gain  error  actuated  flight  control 
systems . 

Professor  Porter's  method  shows  promise  in  enabling  the 
design  of  exceptionally  good  performing  systems  without  many 
of  the  drawbacks  that  occur  in  other  modern  control  methods. 
This  method  leads  to  decoupling  of  the  outputs  and  therefore, 
decoupled  maneuvers  are  readily  achievable.  The  theory 
behind  Professor  Porter's  work  is  presented  in  Chapter  III 
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and  can  also  be  found  in  Reference  3. 

The  longitudinal  and  lateral-directional  dynamics  are 
frequently  considered  separately  in  order  to  reduce  the 
size  of  the  design  problem.  That  procedure  is  used  in  this 
thesis. 

Problem 

Much  of  the  work  on  fighter-type  aircraft  deals  with 
designing  flight  control  systems  considering  the  aircraft 
as  a  rigid  body  vehicle.  These  studies  examine  the  response 
of  the  aircraft  when  performing  various  maneuvers.  It  is 
now  desired  to  include  the  flexible  characteristics  of  a 
fighter-type  aircraft  and  to  examine  the  performance  compared 
to  the  rigid  body  response  when  performing  constant  accelera¬ 
tion,  and  pitch-pointing  maneuvers. 

The  evaluation  of  the  flight  controller  is  defined  in 
this  thesis  as  the  way  the  flexible-fighter  aircraft  deviates 
from  the  rigid  body  response.  The  constant  acceleration 
maneuver  allows  the  pilot  to  command  a  constant  "g"  pull-up. 
The  pitch-pointing  maneuver  allows  the  pilot  to  command  an 
angle  of  attack  without  developing  a  normal  acceleration. 

What  this  means  is  that  the  pilot  can  point  the  nose  of  the 
aircraft  up  or  down  without  changing  its  flight  path. 


Assumptions 

In  order  to  simplify  the  computations  and  analysis  the 
following  assumptions  are  made: 

(1)  The  earth  is  flat  and  non-rotating. 

(2)  The  atmosphere  is  fixed  to  the  earth,  (no  wind) 

(3)  The  acceleration  due  to  gravity  is  constant. 

(4)  The  aircraft  has  a  constant  mass. 

(5)  The  aircraft  is  a  rigid  body  except  for  one  bending 
mode  for  the  longitudinal  axis. 

(6)  The  aircraft  performs  maneuvers  while  flying 
straight  with  wings  level,  unless  otherwise  stated. 

(7)  The  motion  of  the  aircraft  about  a  straight  and 
level  trim  condition  can  be  adequately  modeled  by  a 
linear  dynamical  system  of  equations. 

(8)  The  dynamics  of  the  elevator  and  flaperon  actuators 
can  be  adequately  modeled  by  first-order  differential 
equations. 

These  are  the  basic  assumptions  that  are  made  when 
developing  flight  dynamics  equations  for  a  flexible  fighter- 
type  aircraft  (FFA) . 

Summary  of  Current  Knowledge 

The  aircraft's  aerodynamic  coefficients  of  the  rigid 
body  are  known.  Also,  the  bending  mode  equation  of  the 
flexible  aircraft  is  known.  Pitch-rate  and  normal  accelera¬ 
tion  (Aq)  are  employed  as  command  parameters  for  the  two 
maneuvers  specified  in  the  problem  statement.  When  the 
pilot  demands  a  specific  "g"  level  it  has  a  corresponding 
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pitch  rate  which  can  be  determined  from  the  following 
relationship  (Ref  4): 

q  -  [1845/U]  (An)  (1.3) 

Standards 

In  order  to  have  a  successful  flight  control  design, 
certain  criteria  must  be  established;  therefore,  some  design 
requirements  are  established.  First,  the  design  must 
provide  control  and  stability  augmentation.  Second,  the 
responses  are  to  be  fast  and  well  behaved.  Third,  surface 
deflections  and  rates  are  not  commanded  beyond  preset  limits. 
Fourth,  the  aircraft  rates  and  accelerations  must  be  controll¬ 
able.  Fifth,  the  feedback  signals  can  be  reliably  obtained 
with  sensors.  Last,  the  maneuvers  described  previously  are 
capable  and  available  for  the  flexible  aircraft. 

Approach  Sequence 

The  first  step  is  to  determine  the  bending  mode  equation 
and  to  combine  it  with  the  rigid  body  equations  of  the  fighter. 
This  process  is  described  in  Chapter  II.  Chapter  III  contains 
a  summary  of  the  multivariable  control  design  method  developed 
by  Professor  Brian  Porter  of  the  University  of  Salford,  England. 
Using  this  method,  a  control  design  is  first  developed  for  the 
rigid  body  aircraft.  Then,  the  same  controller  developed  for 
the  rigid  body  is  applied  to  the  flexible-fighter  aircraft 
(FjjkA),  and  the  time  simulation  response  results  are  compared. 
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The  rigid  body  and  the  flexible-fighter  aircraft  responses 
to  the  design  are  examined  in  Chapter  IV.  The  evaluations 
and  conclusions  are  discussed  in  Chapter  V  and  VI , 
respectively. 

There  are  three  appendices  in  this  thesis  that  deal 
with  improvements  which  have  been  incorporated  into  the 
design  program  called  MULTI.  The  changes  given  in  Appendix 
A  are  intended  to  make  the  program  user  friendly  by  making 
it  easy  to  input  corrections  into  the  input  data.  This  is 
accomplished  by  the  Subroutine  Fix.  Appendix  B  helps  the 
user  to  form  combinations  of  states.  For  instance,  to 
look  at  the  flight  path  angle  y  the  equation  6  -  a  must  be 
formed.  This  is  accomplished  by  using  a  new  routine  that 
forms  different  combinations  of  states.  Appendix  C  deals 
with  forming  the  measurment  matrix  M  without  recalculating 
the  Kq  and  K1  matrices  that  may  already  exist.  Also,  this 
improvement  to  MULTI  allows  the  user  to  change  either  or 
both  the  Kq  and  matrices.  The  purpose  of  this  change  is 
to  study  the  possibility  of  achieving  the  desired  tracking 
command  with  a  reduced  number  of  non-zero  elements  in  the 
controller  gain  matrices. 
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II .  Aircraft  Modeling 
Aircraft  Description 

The  vehicle  that  is  described  in  this  chapter  is  a 
flexible-fighter-type  aircraft.  It  is  basically  modeled 
as  a  modified  AFTI/F-16  rigid  body  aircraft  and  is  expanded 
to  include  a  longitudinal  bending  mode  taken  from  a  flexible 
aircraft.  This  has  been  accomplished  by  matching  the 
flexible  aircrafts  dynamic  pressure  (q),  forward  velocity 
(VT),  and  its  speed  relative  to  the  earth  (V  )  with  that  of 
the  AFTI/F-16.  Table  1  shows  three  different  flight  condi¬ 
tions  of  the  AFTI/F-16  and  the  one  for  the  flexible  aircraft 


Table  1 

Matching  Flight  Characteristics 


AIRCRAFT 

FLEXIBLE 

AFTI/F-16 

MACH /ALT. 

0.60/21K 

0.6/20K 

0.8/10K 

0.8/30K 

Po 

0.001225 

0.001267 

0.001267 

0.00089 

vT 

712,6:3 

622.15 

829.54 

795.88 

Ve 

303.1 

269.15 

358.8 

288.66 

q 

311.03 

245.27 

436.03 

282.09 

After  examining  Table  1,  the  flight  condition  that 
most  closely  matches  the  flexible  aircraft  occurs  at  0.8 


mach  at  30,000  feet.  The  next  sections  in  this  chapter 
deal  with  forming  of  the  state  model  for  the  rigid  body 
AFT1/F-16  and  the  steps  taken  to  combine  the  flexible 
aircraft  with  that  of  the  AFTI/F-16.  Reference  4  is  the 
primary  source  used  in  the  next  section. 


Aerodynamic  Model 

For  the  purpose  of  designing  a  digital  flight 
controller,  the  equations  of  motion  must  be  written  as 
state  equations  in  the  form  x  ■  Ax  +  Bu.  The  rigid  body 
space  state  aircraft  model  has  the  form 


— i  m 
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q 

6e 

6f 
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2i 
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X8 

■ 

Xe  -  -gcos0T 

(2.2) 

XA 

= 

X  «  is  c 
u  mU  xu 

(2.3) 

X’ 

. 

X  -  C 

a  m  x 

a 

(2.4) 

a 

xi 

- 

V\  ’  Hr  %  -Bot 

(2.5) 

X' 

«e 

- 
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e  6e 

(2.6) 

®f 

- 

Xfi  "  S2,  cx 

®f  m 

(2.7) 
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(2.1 


(2.1i 

(2.1 


(2.1; 


(2.1! 


Note:  (  )b  denotes  that  the  coefficients  are 

expressed  in  terms  of  the  body  axis  reference 
frame. 

Most  aerodynamic  data  taken  from  wind  tunnel  and  flight 
tests  by  convention,  are  not  in  the  same  reference  frame  or 
in  the  same  units.  For  instance,  by  convention,  static 
derivatives  (C  ,  CT  ,  etc.)  are  listed  in  units  of  per 


degree  (1/deg)  while  the  dynamic  derivatives  (CT  ,CT  etc.) 

q  a 

are  listed  in  units  of  per  radian.  It  should  also  be  men¬ 
tioned  that  the  aerodynamic  derivatives  are  in  the  stability 
axis  reference  frame  while  the  mass  properties  (moments  and 
products  of  inertia,  and  positions  of  sensors)  are  referenced 
in  terms  of  the  body  axis.  Figure  1  illustrates  the  differences 
in  the  body  and  stability  axes.  Therefore,  in  order  for  an 
analysis  to  be  accomplished  the  data  must  be  converted  to  a 
consistent  reference  frame  (Ref.  4). 

The  body  axis  reference  frame  is  chosen  for  this  design. 

The  two  main  reasons  for  choosing  the  body  axis  are  as 
follows : 

1.  The  sensors  providing  feedback  signals  are  actually 
measuring  body  accelerations  and  rates; 

2.  The  pilot  desires  to  control  variables  that  he 
actually  feels,  which  are  body  accelerations  and  rates. 

Because  the  pilot  desires  to  control  body  accelerations 

and  rates,  the  output  equation  takes  the  form 


Z 


u 


0 


at 

»  m 

m 

— 

“  * 

Z  Z 

0 

Z, 

Z. 

6 

a  q 

u 

+ 

5e 

5f 

e 

0  1 

a 

0 

0 

6. 

A 

L 

(2.20) 


The  design  method  used  in  this  report  does  not  provide 
for  a  D  matrix.  Therefore,  the  aircraft  model  must  be 
reformulated  so  that  there  is  no  D  matrix.  This  is  done  by 
augmenting  the  aircraft  state  vector  with  the  actuator  states. 
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Figure  1.  Stability  and  Body  Axis  Differences 


Using  a  first-order  actuator  transfer  function  model  of 
20/(s+20)  for  both  the  elevator  and  flaperon  surfaces  of 
the  AFTI/F-16  results  in  the  following  state  and  output 
equations  that  do  not  contain  a  D  matrix: 


u 

u 
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0 

0 

'  5e 

20 

0 

0 

20 

'CMD 


lCMD 


[0  u  a  q  ie  6fJ 


(2.21) 

T 

(2.22) 


This  form,  where  B  is  in  the  "zero-B2"  form  is  also 

helpful  in  the  design.  One  use  of  this  format  is  that  it 

can  be  used  for  the  direct  evaluation  of  the  transmission 

zeros.  A  relatively  simple  means  of  calculating  transmission 

zeros  is  available  for  systems  in  the  "zero-B2"  form  (Ref.  3). 

The  output  Equation  (2.22)  A„  is  the  longitudinal 

Zcg 

acceleration  in  the  z-direction  at  the  center  of  gravity. 

A  desirable  output  is  the  acceleration  at  the  pilot  station 

(Aq  ).  This  is  accomplished  by  the  following  steps:  (Ref.  4) 
nP 
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Note:  A. 


A  =  A 

D  Q 

P  Cg 


(lx)  (Pr-g)  -  (ly)  (rq+q) 

+  <lz)  (P2+q2) 


(2.23) 


(2.24) 


where  lx,  ly,  and  1  are  distances  in  the  body  from  the 

z 

pilot  station  to  the  center  of  gravity. 

Since  1  and  1  are  much  smaller  than  1  for  the  AFTI/ 
y  z  x 

F-16  aircraft  and  since  the  assumption  is  made  that  the 
longitudinal  and  lateral-directional  motions  can  be  decoupled, 


Equation  (2.24)  reduces  to  the  approximate  form 

V,  *  V. +  (1*)q 

P  eg 

where  lx  ■  13.95. 


(2.25) 


Combining  q  as  given  in  Equation  (2.21),  and  1  with 

A  yields  the  following  equation 
eg 

An  *  8  +  (lx«i-Zu)  U  +  (1vm;-Z„)  a 


x  a  a' 


+  <  1xVZq)  q  +  (1x“je-Z«e)  5e 


*  <1x“Sf-Z{f)  5f 

Therefore,  the  new  output  equation  has  the  form 


(2.26) 


<W  <W-Zo)  <V;*V  (Vq’V  (Vi. -*,  > 


(1xMif-z5,>]  [6  u  “  q  «e  SJT 


(2.27) 
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The  Flexible-Fighter  Aircraft  Model 


The  first  step  in  forming  the  flexible  aircraft  model 
is  calculating  the  A  and  B  matrix  in  Equation  (2.21)  for 
the  rigid  body  aircraft.  The  aerodynamic  data  for  the 
AFTI/F-16  found  in  Table  2  was  developed  by  General  Dynamic 
Corporation  of  Fort  Worth,  Texas  (Ref.  5).  Instead  of 
taking  this  data  and  calculating  the  matrices  by  hand,  a 
computer  program  called  Conversion  and  Transformation  (CAT) 
is  used  (Ref.  3).  The  data  from  the  program  CAT  shown  in 
Table  3  is  put  into  the  state  and  output  Equations  (2.21  & 
2.27) 


• 

0 

0 

0 

0 

• 

U 

U 

0 

0 

• 

- 

a 

=»  A 

a 

+ 

0 

0 

6 

• 

eCMD 

q 

q 

0 

0 

• 

6e 

6e 

20 

0 

^f 

-  CMD. 

_6f 

0 

20 

L 

- 

where 

r  o 

0 

0 

i 

0 

A  ■ 


-33. 15  -0.005647  35.75 

-0.003181  -0.00004481  -0.6488 

0.0003948  0.0005393  3.084 


-43.38 

0.9985 


3.131 

-0.08634 


-0.4555  -11.71 


0 

-1.131 

-0.151 

-3.703 


Table  2 


Aircraft  Data  for  the  0.8  Mach  and  30,000  Feet 
Flight  Condition 


Aircraft  Parameters 

—  2 
q  (dynamic  pressure  -  lbs/ft  )  =  282.33 

2 

S  (wing  reference  area  -  ft  )  =  300.0 
c  (wing  mean  aerodynamic  cord  -  ft)  *  11.32 
b  (wing  span  -  ft)  =  30.0 
VT(trim  velocity  -  ft/sec)  «  795.88 
W  (weight  -  lbs)  =  21,018.0 

Inertias:  I  (slug  -  ft* 2)  =  10,033.4 

I  (slug  -  ft2)  -  53,876.3 

I  (slug  -  ft2)  -  61,278.5 

I„„  (slug  -  ft2)  -  282.132 


Aircraft  Aerodynamic  Coefficients  in  the 
Stability  Axis  System 

a  (deg)  -  3.0615 


V  ‘afe’ 


0.24704 

0.09014 

0.00946 


m 


Cr 


(-As) 


0.01613 


cma  «afe) 

Cm.  (3ig> 

6  1 
(«=) 


0.014623 

0.001886 

-0.011496 


-0.00267 


V3eg> 

c„  (a=r) 


'D. {degJ 

°e 


(a=?> 


0.02746 

0.005118 

0.000084 

0.001015 


Table  3 

0.8  Mach  20.000  Feet  Flight  Condition 


Longitudinal  Body  Axis 
Dimensional  Derivatives 


Z  =  -21018.0 

M  - 

14020.4 

X  - 

-1205.28 

ZQ  -  -675.528 

Ma  * 

1.94806 

Xa  * 

29.753 

Z*  =  -70.228 

6e 

V 

-11.722 

V 

3.1306 

Z,  *  -120.153 

6e 

M.  = 
Si 

-2.7224 

V 

-1.13077 

Zq  *  -2.7683 

Mq  - 

-0.3195 

Xq  = 

0.14806 

Za  «  1.20961 

Ma  = 

-0.13644 

Xa  “ 

-0.064694 

Zu  *  -0.035664 

Mu  “ 

0.0005231 

Xu  * 

-0.006647 

Primed  Dimensional 

Derivatives 

Z^  =  -0.8488 

M*  « 

2.064 

X<x  * 

29.75 

Z6e=  “°* 08824 

M<Se“ 

-11.71 

X6e* 

3.131 

Z£f=  -0.151 

M6f“ 

-2.702 

’  X6f“ 

-1.131 

Z^  *  0.9965 

Mq- 

-.4555 

Xq  = 

-42.38 

Z^  *  -0.00004481 

MU  * 

0.0005292 

XU  * 

-0.006647 

Z£  «  -0.002161 

Me  * 

0.0002948 

xe  = 

-32.15 
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"a 

np 

s 

"0.004113 

0.04305 

704.3 

-3 .585 

-93.12 

82.5“ 

_  q  . 

0 

m 

0 

0 

1 

0 

0 

mi 

X  [e  u  a  q  6e  5{J  T  (2.29) 

The  state  vector  has  the  following  units 


0  (rad.) 
u  (ft/s) 
a  (rad.) 
q  (rad/sec) 
<Sg  (rad.) 

6f  (rad.) 


The  state  space  model  of  a  flexible  aircraft  was  obtained 
from  Dr.  Robert  Schwanz  of  the  Flight  Dynamics  Laboratory  at 
Wright-Patterson  AFB,  Ohio.  The  following  state  quation  is  a 
short  period  approximation  of  the  flexible  aircraft. 


*  *■ 

™  * 

• 

e 

0 

0 

• 

w 

w 

-349.87 

• 

q 

*  A 

q 

+ 

-0.4937 

*i 

*i 

0 

x2 

m  m 

A 

-7328 

] 


where 


0 

0 

1 

0 

0 

A 

-14.3639 

•0.0006139 

-1.09346 

-0.000060164 

8436.64 

-1.6866 

-0.708394 

-0.0001167 

-0.0343707 

-0.0007194 

. 

0 

0 

0 

0 

1 

.16. 76 

-13.96 

7373 

-47.41 

-3.469 

(2.30) 


where  Xj_.  and  X2  represent  the  bending  mode  state  variables. 

Before  a  comparison  can  be  made  between  the  rigid  body 
of  the  AFTI/F-16  with  that  of  the  flexible  aircraft,  w  must 
be  changed  to  a.  This  is  done  by  using  the  equation  a  *  w/u 
where  u  *  713.068  ft/s  and  w  is  in  inches  per  second.  When 
w  is  changed  to  ct  and  the  aircraft  model  is  augmented  with 
the  actuator  dynamics,  the  following  space  state  model  is 
obtained : 


where 


0 

0 

1 

0 

0 

0 

0.001679 

-1.08346 

.9846 

-0.04088 

-0.000002836 

-0.00008243 

-0.0008129 

-4.294 

-1.6588 

-0.49375 

-0.0001168 

-0.00007194 

0 

0 

0 

-20 

0 

0 

0 

0 

0 

0 

1 

1 

16.76 

119400 

7272 

-7328 

-47.41 

-2.469 

(2.31) 

•  • 

Comparing  the  a  and  q  equations  for  both  the  rigid  body 
aircraft  and  the  flexible  aircraft,  the  following  table  is 


formed : 


Table  4 


r 

\ 


Comparison  Between  the  F-16  &  Flexible  Aircraft 


AFIT/F-16  (Rigid  Body) 

Flexible  A/C 

• 

a 

• 

q 

• 

a 

• 

q 

-0.002161 

0.0002948 

0.001679 

-0.0008129 

-0.8488 

2.064 

-1.08346 

-4.294 

0.9965 

-0.4555 

.9846 

-1.6585 

-0.08824 

-11.71 

-0.04088 

-0.49375 

The  magnitudes  of  the  equations  are  nearly  equal  to 
each  other.  As  far  as  modeling  goes  they  are  close  enough, 
because  only  the  bending  mode  equations  of  the  flexible 
aircraft  are  used. 

Combining  the  rigid  body  of  the  AFTI/F-16  and  the 
bending  mode  of  the  flexible  aircraft.  Equation  (2.32)  is 
formed.  The  term  that  is  missing  from  Equation  (2.32)  must 
be  approximated.  To  do  this,  the  following  ratio  is  used: 


0.08824 

-6.151 


7328 

* 


(2.33 


This  yields  *  ■  -12540,  where  represents  the  z-force  with 
either  the  elevator  or  flaperon. 

To  determine  if  the  flexible  aircraft  in  Equation  (2.32) 
represents  a  flexible  fighter-type  vehicle,  the  following 
steps  are  taken: 


(rad) 


Note-  ^represents  the  term  that  is  missing  from  the  flexible  aircraft. 


MX  +  Bx  +  Kx  -  F  (0,a,q,6e,&  6f) 


(2.34) 


where  Equation  (2.34)  represents  a  standard  second-order 
bending  for  the  mass,  spring,  and  damper  system.  Dividing  by 
M,  Equation  (3.34)  becomes 

V  +  I  X  +  |x  -  |  (0,a,q,6e,&6f)  (2.35) 

*X  +  2Cwnx  +  rn^x  *  F' (0,a,q,5e,&6f)  (2.36) 

where  C  is  the  damping  ratio,  u>Q  is  the  undamped  natural 
frequency  and  F'(0,a,q,6e  &  5f)  =  16.76  0  -119400  a  +  7272  q 
-7328  -12540  6f  as  incorporated  in  Equation  (2.32).  For  a 

fighter-type  aircraft  typical  values  for  the  first  bending 
mode  are  5  *  0.025  and  u>n  =  21.0.  Substituting  these  values 
for  c  and  u>n  of  the  fighter  into  Equation  (3.36)  yields  the 


following  equations: 

X  +  1-05  x  +  441  x  -  F,(0,o,q,6e  &  6f)  (2.37) 

Let  Xl  -  X  (2.38) 

x2  "  X1  *  x  (2.39) 

i2  =  -1.05  x2  "441X!  +  F*(0,o,q,«e  &  «f)  (2.40) 


.  .  • 

Substituting  Equation  (2.40)  into  the  place  of  the  x2 
equation  found  in  (2.32),  Equation  (2.41)  is  formed. 

To  determine  if  Equation  (2.41)  actually  represents  a 


1 


flexible  fighter,  its  bending  mode  is  examined  while  performing 

two  different  maneuvers.  The  maneuvers  performed  are  the  g- 

command  pull-up  and  pitch-pointing  which  are  described  in 

Chapter  I.  Before  the  effects  of  structural  bending  on  the 

aircraft  performance  are  looked  at,  columns  7  and  8  of  the  a 

and  q  equations  are  initially  set  to  zero.  Figures  2  and  3 

represents  the  structural  bending  while  performing  the  two 

maneuvers.  Examining  the  two  figures,  the  aircraft's  bending 

is  larger  than  what  would  normally  be  expected.  Therefore, 

the  bending  due  to  the  elevator  and  flaperon  is  reduced  by  a 

factor  of  10.  Performing  the  two  maneuvers  again,  the 

structural  bending  is  illustrated  in  Figures  4  and  5.  Still, 

the  aircraft  bends  more  than  a  fighter  is  expected  to.  For 

instance,  the  g-command  of  2g's  probably  bends  the  aircraft 

no  more  than  a  V  to  1".  Accordingly,  another  value  in  the 
• 

Xg  equation  must  be  reduced.  Knowing  the  bending  equation  of 
a  KC-135,  its  values  are  nearly  of  the  same  magnitudes  as  the 
flexible  fighter,  except  for  the  bending  due  to  a  (Ref.  6). 
Therefore,  reducing  the  influence  on  bending  mode  due  to  a  by 
a  factor  of  10  and  keeping  and  <5f  at  the  reduced  values, 
the  structural  bending  while  performing  the  maneuvers  resembles 
reasonable  values.  Figures  6  and  7  illustrate  the  bending 
taking  place  as  the  aircraft  performs  the  maneuvers. 

Before  a  comparison  is  made  between  the  rigid  body  and 
the  flexible  fighter  an  output  equation  must  be  computed  for 
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'0.00  0.50  1.00  1.50  2.00  2.  SO  3.00  3.50  4.00 

Tine.  SECONDS 


Figure  6.  Structural  Bending  while  performing  the  Pitch- 
Pointing  Maneuver.  (From  Eq.  (2.41)  with  <5e, 
5f>  and  a  reduced) 


r.oo 1.00  IV 


1.00 4.00 1.00 
TWE.  SECONOS 


i.oo  o.o 


Figure  7.  Structural  Bending  while  performing  the  G- 
Command  Pull  Up  Maneuver.  (From  Eq.  (2.41) 
with  ,  5, ,  and  a  reduced) 


the  flexible  aircraft.  Since  the  output  equation  for  the 

rigid  body  has  A  as  an  output,  the  A  equation  for  the 

P  P 

flexible-fighter  aircraft  must  also  be  developed.  This  is 
accomplished  by  using  equations  (2.25),  (2.26),  and  q  in 
Equation  (2.41).  The  following  equation  is  the  result: 


A 


(1^)6  +  <1„M,',-Z„)u  -  (l„M>Z„)a  ♦ 


X  u  u 


x  a  a- 


<1xMq-Vq  +  (1xM«„-V)Se  + 


+  (lxWXl  *  (VW,X2 


e  e 

1  -2 
<2  K2 


(2.42) 


where  Z  -  Z. 


0.  Therefore, 


An  =  [0.004113  0.04305  704.3  -3.585  -92.12  82.5  -0.001628 


-0.001004] 


0 

U 

a 

q 

5e 

6f 

*1 

*2 


(2.43) 


Designing  the  DFSC  for  the  rigid  body  aircraft  and 
comparing  it  with  the  flexible-fighter  aircraft  (FFA),  the 
units  of  the  state  equation  is  changed  to  the  following: 
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6  (deg) 
u  (ft/s) 
a  (deg) 
q  (deg/s) 

5e  (deg) 

Sf  (deg) 

X1  (in) 

X2  (in/s) 

For  the  output  equation  the  units  are 

\  (8) 

P 

q  (deg) 

Tables  5  and  6  are  the  values  for  the  rigid  body  and 
the  flexible-fighter  aircraft,  respectively,  that  are  used 
in  the  designing  of  the  DCFS.  The  flexible-fighter  aircraft 
has  5  different  bending  mode  set  ups  that  are  explained  in 
Chapter  IV. 


TABLE  5 


Rigid  Body  (R.B. ) 


A  MATRIX 


0 

0 

0 

0 

0 

-.5611 

-.006647 

.5192 

-.7396 

.05464 

-.002161 

-.002568 

-.8488 

.9965 

-.08824 

.0002948 

.03032 

2.064 

-.4555 

-11.71 

0 

0 

0 

0 

-20 

0 

0 

0 

0 

0 

X  VECTOR  B  MATRIX 


0  (deg) 

0  0 

u  (ft/s) 

0  0 

a  (deg) 

0  0 

q  (deg/s) 

0  0 

(deg) 

20  0 

6f  (deg) 

-  . 

0  20 

»  _ 

eCDM 

L  CDM 


[ 


2229E-5 

0 


.001337 

0 


C  MATRIX 

.3817  -.001943 

0  1 


-.0547 

0 


OUTPUT 

A  (g) 

P 

q  (deg/s] 


0 

-.01974 
-.151 
-2 . 702 
0 

-20 


(deg) 

(deg) 


04469 
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Ill .  Multivariable  Control  Law  Development 
Introduction 

Of  all  the  digital  control  design  methods  available 
today,  this  thesis  uses  the  method  developed  by  Professor 
Brian  Porter  of  the  University  of  Salford,  England.  A 
summary  of  the  techniques  presented  in  Reference  3  is  given 
in  this  chapter.  The  design  techniques  described  in  the 
references  are  both  conceptually  and  computationally  simple. 

In  addition  to  the  techniques  utilizing  singular  perturba¬ 
tion  methods,  these  are  also  applicable  to  the  design  of 
both  analogue  and  digital  controllers.  The  digital  controller 
developed  in  this  thesis  produces  a  proportional  plus  integral 
control  signal  that  is  held  piece-wise  constant  for  each 
sampling  period. 

Output  decoupling  is  essential  in  the  design  of  many 
Multiple  Input-Multiple  Output  (MIMO)  control  systems.  The 
system  can  be  represented  by  either  a  continuous-time  plant 
model  or  a  discrete-time  plant  model.  For  the  continuous¬ 
time  representation  the  state  and  output  equations  are: 

x(t)  -  A  x(t)  +  B  u(t)  +  D  d(t )  (3.1) 

Z(t)  -  C  x(t)  (3.2) 


where 


A  =  continuous 


-  time  plant  matrix 
B  =  continuous  -  time  control  matrix 
D  »  continuous  -  time  disturbance  matrix 
C  -  continuous  -  time  output  matrix 

The  discrete-time  state  and  output  equations  of  the  plant 
follow  the  form: 

x(kT+T)  =  <J>x(kT)  *  'i'u(kT)  +  Ad(kT)  (3.3) 

£(kT)  =  Tx(kT)  (3.4) 

where 

<J»  *  sampled-data  plant  matrix 
»  sampled-data  control  matrix 
A  «  sampled-data  disturbance  matrix 
T  «  sampled-data  output  matrix 


and 

<t>  ■  exp  (AT)  (3.5) 

V  «  exp (AT)  B  dt  (3.6) 

A  -  exp (AT)  D  dt  (3.7) 

r  »  C  (3.8) 


The  dimensions  n,  m,  and  p  represent  the  number  of  states, 
inputs,  and  outputs,  respectively. 

The  plant  in  state  space  form  has  the  following  format: 
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A 

-11 

-12 

*1 

Si 

= 

—  — — -f - - 

+ 

—  — 

-22 

A 

-22 

1  - 

r 

IX 

(O 

1 _ 

—2 
m  m 

(3.9) 

(3.10) 


where  B2  is  a  square  matrix  with  a  row  and  column  size  equal 
to  the  number  of  inputs  in  the  vector  u;  and  £2  is  a  square 
matrix  with  the  number  of  rows  and  columns  equal  to  the 
number  of  outputs.  Also,  A22  has  the  same  dimensions  as 
and  £2 .  Knowing  these  dimensions  allows  the  system  matrices 
to  be  easily  partitioned  as  shown  above.  It  is  preferable, 
but  not  necessary,  if  «  0.  This  may  occur  as  in  a  kinema¬ 
tic  equation  or  may  be  achieved  by  a  transformation  of 
variables. 

Requirements  to  have  a  successful  design  using  Porter's 
techniques  are  as  follows: 

(1)  The  system  must  be  controllable  and  observable. 

(2)  The  number  of  inputs  must  equal  the  number  of  out¬ 
puts. 

Two  design  methods  discussed  in  this  chapter  are  the 
regular  and  the  irregular  designs  for  the  case  where  the 
state  and  output  equations  are  known.  The  choice  between 
the  two  design  methods  is  determined  by  the  type  of  plant 
being  considered,  and  the  outputs  which  are  to  be  controlled. 


Regular  Plants 

To  be  considered  as  a  "regular"  plant,  the  rank  of 
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the  first  Markov  parameter,  [  CB  ]  ,  must  have  full  rank. 

Because  [CB]  is  of  full  rank,  the  gain  matrix  for  feedback 
control  is  found  using  the  inverse  of  [CB], 

The  control  law  is  a  proportional  plus  integral  out¬ 
put  feedback  law  expressed  in  the  descrete  case,  shown  in 
Figure  8  (for  the  regular  plant  the  M  matrix  is  not  utilized), 
as : 

u(kT)  -  (jjr)  [Kq£ (kT)  +  K^z(kT)]  (3.11) 

where 

i  -  is  a  scalar  gain  which  is  equal  to  the  sampling 
frequency 

K«  =  is  the  proportional  gain  matrix  for  the  error 
signal  e(kT) 

K-  -  is  the  gain  matrix  for  the  backward  difference 
‘  of  the  error  signal,  which  is  designated  as  z(kT) 

Although  a  continuous  integrator  is  shown  in  Figure  8,  the 

digital  implementation  represents  the  discrete  state  equation 

jg(kT)  =  r(kT)  -  y(kT) 
z  [  (k+l)T  ]  *  z(kT)  +  Te(kT) 

For  the  continuous  case,  shown  in  Figure  9,  the  law  is 
expressed  as 

u(t)  -  g(KQ©(t)  +  K1/e(t))dt  (3.12) 


35 


where 


g  is  the  forward  path  gain 

Kq  is  the  proportional  gain  for  the  error  signal 
is  the  gain  for  the  integral  of  the  error  signal 

The  Kq  and  the  matrices  are  determined  using  the  following 
equations : 

Kq  =  (CB)-1Z  (3.13) 

Kx  =  a(CB)-1Z  (3.14) 

where 

a  is  the  ratio  of  the  error  e  to  the  integral  of 
error  z 

_!  is  the  diagonal  weighting  matrix  whose  scalar  elements, 
cri ,  are  chosen  by  the  designer .  These  values  determine 
the  respective  gains  or  weighting  of  each  error  signal 
on  the  control  surface. 

The  error  vector  e  is  expressed  as  e  «  v-jr.  Each  command 
in  the  v  vector  is  compared  with  its  corresponding  output 
element  in  the  £  vector  and  is  the  input  for  the  proportional 
plus  integral  controller.  Since  the  complete  system  includes 
a  vector  integrator,  the  steady-state  value  of  the  error 
vector  is  zero  for  a  constant  command  vector  input,  and 
tracking  is  achieved.  Therefore,  because  the  error  vector  is 
driven  to  zero,  the  resulting  output  equals  the  command  input. 

Irregular  Plant 

Plants  in  which  the  First  Markov  parameter,  [  CB],  is 
rank  deficient  are  designated  as  "irregular".  For  the 


I 


i 

t 


I 
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irregular  plant  the  M  matrix  in  Figure  8  is  utilized  to 
allow  the  design  to  be  accomplished.  An  F  matrix  is  formed, 
as  shown,  and  used  in  place  of  the  C  matrix: 


I  =  [F±  j  F2J  (3.15) 

where 

II  =  l£x  +  MxiJ  (3 . 3.6) 

F2  -  [£2  +  MA12]  (3.17) 


The  elements  of  the  measurement  matrix  M  are  selected  to 
obtain  a  matrix  £FB]  having  full  rank  and  thus  being  invert¬ 
ible.  M  is  sleeted  to  be  as  sparse  as  possible  (Ref.  7). 

The  control  law  for  the  discrete  case  in  Equation  (3.11) 
and  for  the  continuous  case  in  Equation  (3.12)  are  the  same 
as  the  regular  plant,  except  that  the  Kq  and  the  matrices 
are  determined  by  the  following  equations: 


Ko  *  (FB)"1! 


(3.18) 


K.i  *  <X(£B)-1I 


(3.19) 


In  the  irregular  design  the  e  vector  is  defined  by 
e  ■  v-w  (3.20) 

where 

w  -  jr  +  Mx^  (3.21) 
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These  outputs  w  are  developed  with  the  measurement  matrix 
M  and  the  derivatives  x^  of  the  system  states.  The  error 
vector,  e,  for  the  irregular  design  becomes  the  same  as 
the  regular  design  when  the  states  reach  their  steady-state 
value.  This  can  be  seen  as  follows:  With  a  constant 
command  vector  v  the  states  reach  constant  steady-state 
values  and  therefore: 

• 

Xj  =  0  (when  x^  reaches  steady-state) 
so  that 

^ss  "  Zss  (3’22) 

Substituting  Equation  (3.22)  into  Equation  (3.20)  the 
steady-state  error  vector  again  becomes  e  «  v-£.  Therefore, 
since  the  error  e  *  0,  as  before,  the  output  equals  the 
command  input,  that  is  £  ■  v. 

Asymptotic  Characteristics 

The  system  transfer  function  G(X)  »  £(^I.-A)“1B  assumes 
the  asymptotic  form 

£(X)  -  £(X)  +  £(X)  3.23) 

where 

£(X)  is  the  "slow"  transfer  function  matrix, 

A 

£(X)  is  the  "fast"  transfer  function  matrix. 


and 


The  system  reaches  these  asymptotic  values  as  the  gain 
factor  g->-»  for  the  continuous  design  or  as  the  gain  factor 
^►09  for  the  discrete  case.  The  poles  of  £(X)  consist  of 
two  sets  of  roots  that  are  designated  as  and  Z2<  The 

A 

poles  of  the  fast  transfer  function  T(X)  are  designated  as 
Zg.  For  a  regular  design  the  slow  modes  associated  with 
the  set  of  roots  Z. ^  become  uncontrollable  while  modes 
associated  with  the  set  of  roots  Z 2  (composed  of  the  trans¬ 
mission  zeroes)  become  unobservable  as  the  gain  factor 
increases  toward  infinity.  As  this  happens,  the  system 
response  is  increasingly  dominated  by  the  fast  modes  which 
are  associated  with  the  set  of  roots  Zg.  Thus,  the  slow 
modes  contained  in  £(X)  disappear  from  the  overall  transfer 

A 

function  £(X),  and  only  the  fast  modes  in  £(X)  remain  for  a 
regular  design.  The  expressions  for  the  asymptotic  modes 
are  given  in  Table  7. 

The  system  transfer  function  G(X)  approaches  a  diagonal 
asymptotic  form  of  increasingly  non-interacting  control. 
Therefore,  as  a  result,  the  output  responses  become  decoupled 
as  the  gain  increases.  The  asymptotic  closed-loop  transfer 
function  of  the  discrete  system  has  the  form: 

I(X)  -  dlag  . (3.24) 

For  the  continuous  system,  the  asymptotic  transfer  function 
has  the  form 
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g0l  ga2  gai 

r_(X)  -  diag  X+goT'  X+go„' - '  X+ga .  (3.25) 

where  the  a's  come  from  the  weighting  matrix  Z. 

For  the  irregular  plant,  the  previous  equations  can  be 
applied.  The  equations  in  Table  5  are  the  same  for  the 
irregular  design  after  replacing  the  C  matrix  with  the  F 
matrix.  The  component  is  replaced  by  F^  and  is 
replaced  by  Fg  as  calculated  in  Equations  (3.16)  and  (3.17). 

It  should  be  noted  for  the  irregular  case  the  CQ  is  different 
as  indicated  in  the  table.  Also,  for  the  irregular  design, 
as  the  slow  modes  associated  with  the  set  of  roots  become 
uncontrollable  the  modes  associated  with  the  set  of  roots  Z2 
remain  observable  as  the  gain  factor  increases  toward 
infinity . 

The  asymptotic  closed-loop  system  roots  can  be  calculated 
from  the  equations  in  Table  5  and  are  composed  of  the  sets 
&1 ,  Sg,  and  The  following  rules  can  be  used  to  estimate 

root  locations  in  the  s-plane  as  g-*-®. 

The  set  of  roots  Z^  are  m  in  number  and  are  all  equal  to 
5,  the  ratio  of  integral  to  proportional  gain  values,  if  = 
agg.  The  set  of  roots  _g2  are  n-m  in  number  and  located  at  the 
transmission  zero  locations.  Thus,  the  transmission  zeros 
(Zt)  can  be  calculated  from  the  Z2  equation  in  Table  7  if 
the  B  matrix  has  the  "zero-B2"  form  shown. 

The  set  of  roots  are  m  in  number  and  for  the  continuous 
case  they  are  located  at  -ga^  on  the  real  axis  where  g  is  the 
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IV.  Responses  for  the  Rigid  Body  and  the 


Flexible-Fighter  Aircraft 
Frequency  Responses 

The  effect  of  a  bending  mode  on  the  frequency  response 

of  the  longitudinal  acceleration  at  the  pilot's  station 

(Afl  )  with  either  an  elevator  or  flaperon  input  is  examined 

in  this  section.  Figures  10  and  11  show  the  frequency 

responses  A  /6  and  A  /6fl  respectively,  for  the  rigid 
P  P 

body  aircraft.  Figures  12  and  13,  correspond  to  the 
flexible-fighter  aircraft  (FFA)  as  modeled  in  Table  6, 
Chapter  II. 

A  series  of  five  representations  of  the  bending  mode 

and  its  coupling  into  the  aircraft  motion  variables  is 

studied.  This  permits  a  study  of  the  effects  of  a  wide 

range  of  aircraft  structural  bending  characteristics  on 

the  aircraft  tracking  performance.  Examining  the  effect 

that  the  bending  mode  has  an  A  ,  the  factors  dealing  with 

P 

the  bending  mode  states  and  X2  are  increased  by  a 
factor  of  10  which  yields  the  following  output  equation: 

A  »  [0.000002229  0.001337  0.3817  -0.001943  -0.0547 

p 

0.04469  -0.0005056  -0.0003118]  x 

(4-1) 

where  x  ■  [0  u  a  q  5f  x2]T 

Figures  14  and  15  show  the  frequency  responses  of  A  /S 

p 


FIGURE  10.  RN/OELTR  ELV.  FOR  THU  RIGIO  800r  ( R.B.) 


FIGURE  11.  RN/OELTfl  FPN.  FOR  THE  R.B 


MOttlTUOC  I  DECIBELS  >  |  nftOUlTUOE  IOCCIKLS) 


FIGURE  14.  RN/DELTfl  ELV.  FOR  THE  FFfilG. 
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and  A  /$.,  respectively,  using  Equation  (4-1),  which  are 
P 

represented  as  FFA10 . 

To  further  examine  these  effects,  the  and  X2 
coefficients  are  increased  by  a  factor  of  100  which  is 
shown  in  Equation  (4-2): 


Aq  *  [0.000002229  0.001337  0.3817  -0.001943  -0.0547 

P  0.04469  -0.005056  -0.003118]  x 


Figures  16  and  17  are  the  results  using  the  output  Equation 
(4-2).  The  new  aircraft  model  is  represented  as  FFA100. 

The  output  Equation  A  is  originally  computed  by 
Equation  (2.25)  as  A  +  1  (q) ,  as  described  in  Chapter  II. 

Cg 

As  a  further  modification  of  the  aircraft  model  the  and 

X2  factors  found  in  the  q  equation  are  increased  by  the 

same  amount  as  in  the  output  equation.  Figures  18  and  19, 

represented  as  FFAlOq,  are  the  responses  AQ  /6fi  and  Aq  /6f, 

P  P 

respectively,  using  Equation  (4-1)  for  A  and  increasing 

P 

the  x^  and  X2  factors  in  the  q  equation  by  10.  Figures  20 

and  21,  are  the  results  of  increasing  the  and  X2  factors 

•  •  ■ 

in  the  q  equation  by  100  and  using  the  output  Equation  (4-2) 

for  A  .  These  figures  represent  the  aircraft  model  FFAlOOq. 
P 

Time  Responses 

This  section  examines  the  digital  f'u^f.c  control  design 


for  the  rigid  body  and  for  the  five  different  flexible- 
fighter  models  performing  the  q-command  pull-up  maneuver 


and  the  pitch-pointing  maneuver.  The  five  flexible- 
fighter  models  are  those  described  in  the  first  sections 
of  this  chapter. 

The  following  design  data  is  used  to  develop  the 
simulation  time  responses: 


Design  Data  for  the  G-Command  System 
Flight  Condition:  0.8  Mach  30,000  feet 

T  =  0.02 

Computational  Delay:  0 
a  *  0.2 

To.  025  0 

0  2.0 
2.0 


Rigid  Body  Model 

=  f°  0 

lo  0 

1-0.002013  -0.02664*1 

,008724  -0.03260J 

1-0.01007  -0.1332 

.04362  -0.1630 


s  -  r°  °  °  °  i 

|o  o  o  o.ij 

Ko  -  P0J 

0  Lo.< 

Ki  *  P0,1 

1  Lo.< 


r 


Flexible-Fighter  Models 

0  0  0  0  0  o“I 

0  0  0  0.1  0  oj 


] 


Input  ramp-up  time:  0.4  seconds 

Input  commands:  A  -  lg  step 

P 


_  _  1845  » 

<!  -  “C  An 


2.32  deg/sec 


P 
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NOTE:  Step  inputs  are  "ramped-up"  to  the  steady- 
state  value  over  a  specified  time.  This 
time  is  designated  as  the  "input  ramp-up 
time"  given  above. 


Design  Data  for  the  Pitch-Pointing  System 
Fight  Condition:  0.8  Mach  30,000  feet 


T  =  0.02 

Computational  Delay :  0 

a  =  0.5 

I  =  T1*5 

L°  1*°. 

e  =  2.5 


Rigid  Body  Model 
M 


Flexible-Fighter  Models 


0  0  0  0  ' 

0  0  0  0  0  o“ 

* 

M  - 

0  0  0  0.1 

0  0  0  0.1  0  0_ 

■  [t 
*  ■  [: 


-0.3774 


636 

-0.7549 
272 


-0.04162' 
-0 . 05094_ 


-0.08324 

-0.1019 


'] 


Input  ramp-up  time:  0.0  seconds 
Input  command:  AQ  ■  0.0 

q  -  1.0  deg/sec  step 
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NOTE:  Step  commands  are  ramped-up  to  the  steady- 

state  value  over  a  specified  time.  This  time 
is  designated  as  the  "input  ramp-up  time" 
given  above. 

NOTE:  The  pitch-pointing  maneuver  is  usually  thought 
of  as  a  pilot  commanding  a  certain  angle,  0, 
without  changing  its  flight  path.  Actually, 
the  pilot  commands  a  pitch  rate,  q,  and  when 
it  reaches  the  angle,  0,  the  pitch  rate  is 
brought  to  zero.  This  is  the  reason  that  the 
design  for  the  pitch  rate  maneuver  commands  a 
pitch  rate,  q,  and  not  an  angle,  6. 

Figures  30a-d  and  31a-e  are  the  results  of  commanding  one 

degree  of  pitch-pointing  for  the  rigid  body  and  the  flexible- 

fighter  aircraft  (FFA) .  That  is: 


Input  Ramp-up  time:  0.0  seconds 

Input  command:  A  *  0.0 

P 

q  “  1.0  deg/sec  pulse  for  1  sec  duration 
Maximum  Command  Responses 

This  section  contains  responses  to  maximum  estimated 
commands  which  have  been  used  to  select  system  gains  so  that 
control  surface  position  limits  are  not  exceeded.  The  g- 
command  system  design  includes  only  the  rigid  body  and  the 
flexible-fighter  model  as  stated  in  Table  6  of  Chapter  II. 


Figure  22a.  G-Command  System  Response  to  a  lg  Step 

Command  for  the  Rigid  Body  Aircraft.  (0.8 
Mach  30,000  ft) 


Figure  22b.  G-Command  System  Response  to  a  lg  Step 

Command  for  the  Rigid  Body  Aircraft.  (0.8 
Mach  30,000  ft) 
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Figure  22c.  G-Command  System  Response  to  a  lg  Step 

Command  for  the  Rigid  Body  Aircraft.  (0.8 
Mach  30,000  ft) 


Figure  22d.  G-Command  System  Response  to  a  lg  Step 

Command  for  the  Rigid  Body  Aircraft.  (0.8 
Mach  30,000  ft) 
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Figure  24a.  G-Command  System  Response  to  a  lg  Step 
Command  for  the  FFA,  for  the  FFA10,  and 
for  the  FFAlOq.  (0.8  Mach  30,000  feet) 


Figure  24b.  G-Command  System  Response  to  a  lg  Step 
Command  for  the  FFA,  for  the  FFA10,  and 
for  the  FFAlOq.  (0.8  Mach  30,000  feet) 


Figure  24c.  G-Command  System  Response  to  a  lg  Step 
Command  for  the  FFA,  for  the  FFA10,  and 
for  the  FFAlOq.  (0.8  Mach  30,000  feet) 


Figure  24d.  G-Command  System  Response  to  a  lg  Step 
Command  for  the  FFA,  for  the  FFA10,  and 
for  the  FFA  lOq.  (0.8  Mach  30,000  feet) 


Figure  24e.  G-Command  System  Response  to  a  lg  Step 
Command  for  the  FFA,  for  the  FFA10,  and 
for  the  FFA  lOq.  (0.8  Mach  30,000  feet) 


Figure  25a.  Pitch-Pointing  System  Response  to  a  1  deg/sec 
Step  Command  for  the  FFA ,  for  the  FFA10  and 
for  the  FFAlOq.  (0.8  Mach  30,000  feet) 


Figure  25b.  Pitch-Pointing  System  Response  to  a  1  deg/sec 
Step  Command,  for  the  FFA,  for  the  FFA10,  and 
for  the  FFAlOq.  (0.8  Mach  30,000  feet) 


1.50  2.00  2.50 

TINE.  SEC0N0S 


3.50  4.00 


Figure  27a.  Pitch-Pointing  System  Response  to  a  1  deg/sec 
Step  Command  for  the  FFA100.  (0.8  Mach 
30,000  feet) 


t.H  2.00  2.10 

HUE.  SECONDS 


Figure  27b. 


Pitch-Poiating  System  Response  to  a  ideg/sec 
Step  Command  for  the  FFA100.  (0.8  Mach  30,000 
feet) 
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| _ _ TIME.  SECONDS _ | 

Figure  27c.  Pitch-Pointing  System  Response  to  a  1  deg/sec 
Step  Command  for  the  FFA100.  (0.8  Mach  30,000 
feet ) 


Voo  <Tio  Too  ITio  Too  Tio  Too  Tso  *■  oo 
_ TINE.  SECONDS _ _ 

Figure  27d.  Pitch-Pointing  System  Response  to  a  1  deg/sec 
Step  Command  for  the  FFA100.  (0.8  Mach  30,000 
eet ) 


Figure  28a.  G-Command  System  Response  to  a  lg  Step 

Command  for  the  FFAlOOq.  (0.8  Mach  30,000 
feet) 


L*tf.oo  i'.oo  Too  Too  Too  Tod  Too  Too  Too 

_ HUE.  3EC0N0S _ 

Figure  28b.  G-Command  System  Response  to  a  lg  Step 

Command  for  the  FFAlOOq.  (0.8  Mach  30,000 


Figure  28c.  G-Command  System  Response  to  a  lg  Step 

Command  for  the  FFAlOOq .  (0.8  Mach  30,000 
feet) 


Figure  28d.  G-Command  System  Response  to  a  lg  Step 

Command  for  the  FFAlOOq.  (0.8  Mach  30,000 
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Figure  30c.  Pitch-Pointing  Command  System  Response  to  a 
1  deg/sec  pulse  of  1  sec  duration  for  the 
Rigid  Body  Aircraft. 


Figure  30d.  Pitch-Pointing  Command  System  Response  to  a 
1  deg/sec  pulse  of  1  sec  duration  for  the 
Rigid  Body  Aircraft. 


Figure  31e.  Pitch-Pointing  Command  System  Response  to  a 

1  deg/sec  pulse  of  1  sec  duration  for  the  FFA. 


The  responses,  Figures  32a-d  and  33a-e,  of  this 
section  are  not  intended  to  indicate  adequacy  of  both 
the  rigid  body  or  the  flexible-fighter  models  in  properly 
representing  such  large  maneuvers.  They  are  simply  used 
as  a  guideline  in  establishing  realistic  gains.  The 
design  parameters  given  in  the  last  section  apply  to  these 
responses  with  the  exception  that  the  inputs  have  been 
appropriately  increased. 

Computational  Delay  Responses 

This  section  contains  the  time  responses.  Figures  34 
through  37,  of  the  design  that  has  a  computational  delay 
incorporated.  The  time  delay  in  the  design  represents  the 
expected  microprocessor  computational  time  delay  in 
generating  the  control  law.  The  simulation  time  responses 
are  for  both  the  rigid  body  and  the  flexible-fighter 
aircraft  (FFA). 

The  design  parameters  are  the  same  as  stated  in  the 
second  section  except  that  there  is  a  computational  delay 
of  l  sampling  period. 

The  evaluation  and  conclusions  that  are  drawn  from  the 
design  are  discussed  in  Chapter  V  and  VI,  respectively.  The 
evaluation  is  based  on  a  pilot's  expected  rating  for  the 
flight  controller. 
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Figure  32c.  G-Command  System  Response  to  a  6.5g  Step 

Command  for  the  Rigid  Body.  (0.8  Mach  30,000 
feet) 


G-Command  System  Response  to  a  6.5g  Step 
Command  for  the  Rigid  Body.  (0.8  Mach  30,000 
feet) 


Figure  32d. 


Figure  33c.  G-Command  System  Response  to  a  6.5g  Step 

Command  for  the  FFA.  (0.8  Mach  30,000  feet) 


Figure  33d.  G-Command  System  Response  to  a  6.5g  Step 

Command  for  the  FFA.  (0.8  Mach  30,000  feet) 
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Figure  34a.  G-Command  System  Response  to  a  lg  Step 

Command  with  a  Computational  Delay  of  1  for 
the  Rigid  Body.  (0.8  Mach  30,000  feet) 


■es  ra •• 
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Figure  34b.  G-Command  System  Response  to  a  lg  Step 

Command  with  a  Computational  Delay  of  1  for 
the  Rigid  Body.  (0.8  Mach  30,00  feet) 
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Figure  34c.  G-Command  System  Response  to  a  lg  Step 

Command  with  a  Computational  Delay  of  1  for 
the  Rigid  Body.  (0.8  Mach  30,000  feet) 


Figure  34d.  G-Command  System  Response  to  a  lg  Step 

Command  with  a  Computational  Delay  of  1  for 
the  Rigid  Body.  (0.8  Mach  30,000  feet) 


Figure  35a.  Pitch-Pointing  Command  System  Response  to 

a  1  deg/sec  Step  Command  with  a  Computational 
delay  of  1  for  the  Rigid  Body.  (0.8  Mach 
30,000  feet) 


Figure  35b.  Pitch-Pointing  Command  System  Response  to  a 
1  deg/sec  Step  Command  with  a  Computational 
delay  of  1  for  the  Rigid  Body.  (0.8  Mach 
30,000  feet) 
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Figure  35c.  Pitch-Pointing  Command  System  Response  to  a 
1  deg/sec  Step  Command  with  a  Computational 
delay  of  1  for  the  Rigid  Body.  (0.8  Mach 
30,000  feet) 
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Figure  36c.  G-Command  System  Response  to  a  lg  Step 

Command  with  a  Computational  delay  of  1  for 
the  FFA.  (0.8  Mach  30,000  feet) 


Figure  36d.  G-Command  System  Response  to  a  lg  Command 

with  a  Computational  delay  of  1  for  the  FFA. 
(0.8  Mach  30,000  feet) 


TIME.  SECONDS 


Figure  36e.  G-Command  System  Response  to  a  lg  Command 

with  a  Computational  delay  of  1  for  the  FFA. 
(0.8  Mach  30,000  feet) 
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'a.  Pitch-Pointing  System  Response  to  a  1  deg/sec 
Step  Command  with  a  Computational  Delay  of  1 
for  the  FFA.  (0.8  Mach  30,000  feet) 


b.  Pitch-Pointing  System  Response  to  a  1  deg/sec 
Step  Command  with  a  Computational  Delay  of  1 
for  the  FFA.  (0.8  Mach  30,000  feet) 
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Figure  37c.  Pitch-Pointing  System  Response  to  a  1  deg/sec 
Step  Command  with  a  Computational  Delay  of  1 
for  the  FFA.  (0.8  Mach  30,000  feet) 


Figure  37d.  Pitch-Pointing  System  Response  to  a  1  deg/sec 
Step  Command  with  a  Computational  Delay  of  1 
for  the  FFA.  (0.8  Mach  30,000  feet) 


V.  Pilot  Rating  Evaluation 

Neal-Smith  Criterion 

After  a  digital  flight  controller  is  designed  a 
means  of  determining  its  acceptability  by  a  pilot  is 
needed.  One  method  of  predicting  a  pilot's  acceptance  for 
a  longitudinal  compensatory  tracking  task  has  been  suggested 
by  Peter  Neal  and  Roger  Smith.  The  complete  theoretical 
basis  and  the  data  used  for  the  development  of  this  criter¬ 
ion  are  provided  in  Reference  8.  A  summary  of  the  Neal- 
Smith  Criterion  for  predicting  pilot  rating  is  presented  in 
this  section. 

An  assumption  is  made  when  using  the  criterion;  that 
is,  the  acceptability  of  an  aircraft's  maneuvering  response 
in  performing  a  specific  task  can  be  expressed  in  terms  of 
a  function  representing  a  pilot's  compensation  needed  to 
achieve  some  "minimum  standard  of  performance"  with  the 
least  tendency  to  develop  a  pilot  induced  oscillation  (FIO) . 
The  standard  of  performance  is  established  by  the  maneuver 
which  the  pilot  is  performing.  For  fighter  aircraft  the 
precise  control  of  pitch  attitude  is  assumed  to  be  the 
critical  task  in  the  reference.  The  control  of  pitch 
attitude  is  modeled  as  shown  in  Figure  38.  The  model  and 
definitions  presented  in  the  figure  are  direct  extracts 
from  Reference  8. 

In  order  to  use  the  Neal-Smith  criterion,  the  engineer 
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Using  this  model,  the  following  terms  ere  defined: 


Bandwidth  (BW);  Bandwidth  is  defined  is  the  frequency  for 

which  the  closed-loop  Bode  phase.  ,  isequaf 

to  -90  degrees.  It  is  a  measure  of  how  quickly  the  pilot 
can  move  the  airplane's  nose  toward  the  target. 


Droop;  Droop  is  defined  as  the  maximum  excursion  of 

closed-loop  Bode  amplitude.  |®/^|  •  below  the  0  dB 


,11ns  for  frequencies  less  than  BW.  In  the  absence  of  large 
oscillations,  droop  is  a  measure  of  how  slowly  the  nose 
settles  down  on  target. 


Standard  of  Performance:  A  minimum  bandwidth.  (BW) 

nun 

of  I.  5  rad/sec,  and  a  maximum  droop  of  3  dB : 


*  ?/*«) 
and  !•/•*! 


greater  th-.n  (-90)  degree 
greater  than  (-3)  dB 


) 


fo.*  CD  less 
than  3. 5 


PIO  Tendency;  The  ten  dency  to  oscillate  cr  PIO  is 

deiined  in  terms  of  the  Bede  magnitude  of  any  closed-loop 

resonant  peak,  |0/dt|  max  »  that  results  from  the  pilot's 

efforts  to  achieve  the  performance  standards. 


Pilot  Compensation;  The  pilot's  physical  and  mental  workload 
required  to  achieve  the  standard  of  performance  is  defined 

in  terms  of  the  phase  of  his  compensation  at  a)  «  (BW)^ 


Figure  38.  Neal-Smith  Model  for  Pitch  Attitude  Control 


selects  values  for  Kp,  Tp1»  and  tp2  in  the  pilot  compensation 
model  that  yield  the  smallest  maximum  value  for  1  © / ©  [  for 
the  closed-loop  system  while  maintaining  a  minimum  bandwidth 
of  3.5  radians  per  second,  and  a  maximum  droop  of  3dB.  The 
procedure  for  selecting  the  values  are  covered  later. 

The  Nichols  chart  is  used  to  determine  the  maximum 
value  of  1 9/6  |  and  the  phase  angle  at  the  frequency  of  3.5 
rad/sec  for  each  of  the  pilot  model  parameter  selections. 

The  value  of  Kp  and  the  pilot  compensation  angle  determined 
by  the  xpl  and  Tp2  values  are  plotted  in  Figure  39.  This 
figure  is  used  to  evaluate  the  design,  that  is,  to  determine 
a  predicted  pilot  rating  (PR)  and  the  acceptability  of  the 
design. 

To  apply  this  method,  the  airframe  plus  flight  control 
system  transfer  function,  9/Fs,  must  first  be  determined 
for  the  longitudinal  system,  where  Fg  is  the  stick  force 
applied  by  the  pilot.  Equation  (5.1)  represents  the  control 
system  transfer  function  matrix  G(s). 


Application 

Since  the  normal  acceleration  A„  and  q  are  related 

n 
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Figure  39.  Neal-Smith  Criterion  for  Fighter  Maneuvering  Dynamics 


quantities,  both  inputs  must  be  considered  in  the  pitch 
attitude  output.  Figure  40  relates  stick  force  input  (Fe) 
to  the  control  commands  v^  and  v2  for  both  the  g-command 
and  pitch-pointing  maneuvers.  For  the  g-command  maneuver 

vx  =  Fs  and  v2  «  (i§^)(Ffl)  (5.2) 

Thus 


2.  (Fe) 
v^  v  s' 


.1845 


)(FC) 


and 


(5.3) 


For  the  pitch-pointing  maneuver 


(5.4) 

(5.5) 


Vj  -  0.0  and  v2  -  Fg 

thus 


(5.6) 


I  (a  > 

s  W2> 


(5.7) 


Equations  (5.8)  and  (5.9)  are  the  transfer  function  0/F_, 
obtained  for  Equation  (5.5),  of  the  rigid  body  and  the 
flexible-fighter  aircraft,  respectively,  performing  the 
q -command  maneuver. 


CO  3 

a| 

£  C 

y 

p> 


o 

o 

i 

o 


■5  E  I 

■p  0)  I 

C  -P  f 

•H  0]  • 
0  >» 
a  co 


x: 

H 

■P 

0) 

•H 

£ 


101 


E 

•p  a> 
o  -p 

rH  (0 
•H  >» 

a  co 


£  be 
•p  c 


as  -p 
e 
E  "H 
<D  O 

■p  a 
« 

>»  JS 

co  o 

•p 
TJ  -H 

c  a 

CS 


tc 

e 


rH  73 
«  3 
C  rH 
•P  O 

■o  c 

3  rH 
•P 

•H  T3 
be  c 
c  ea 
O 

J  c 
o 

<H  -H 
O  -P 


E  73 

a  T3 
u  < 
bo 

0$  <D 
•H  o 
Q  as 


.*  s- 
y  a> 
o  -p 

rH  C 

CQ  rn 


O 


0) 

U 

3 

bfi 


1 


0.  =  927.3(s+0. 00036+10,0034^ (s+0. 0094) . 

Fg  s ( s-0 . 00075 ) ( s+0 . 0037 )  ( s+0 .0073 ) ( s+1 . 012 ) 

(s+0. 986) (s+5.0) _ 

(S+5.79+J2.54) (s+48.82)  (5.8) 

9  _  927. 3(S+0.521+.120, 99)  (S+o.QQ88±,iO. 00874)  . 

Ps  s ( s+0 . 528+ J20799 ) ( s-0 . 00155+ jO7O075 ) (S+6 . 6l354 ) 

(s-0. 0076) (s +0.984) (s*5.0) 

(s+1.01) (S+5.79+j2.54) (S+48.82)  (5.9) 

Equations  (5.10)  and  (5.11)  represent  the  transfer  functions 
of  0/Fg,  obtained  from  Equation  (5.7),  of  the  rigid  body  and 
the  flexible-fighter  aircraft,  respectively,  for  the  pitch¬ 
pointing  maneuver. 


e 
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625.0(S+2.0) (s+l. 818) (S+Q. 0025+jO. 00703) 
s  s ( s+72 . 8 ) ( s+8 . 12 ) ( s+i . 965 ) ( s+1. 898 ) ( s+0 . 0037+ jO . 005 ) 


(S+0. 016) 
(s+0. 0139) 


0  625.0 ( s+0. 5242±J20. 99) (S+2.0)( s+1. 817) ( s-0. 0075) 

F  s( s+0 .529+ j20 . 99 ) ( s+72 .8 ) ( s+8 .12) ( s+1. 966 ) ( s+1 . 8 

s  ™ 

( S+0. 014±j0. 0123) _ . 

( s+0. OSOI+JO. 0234) (S-0. 0196) 


(5.10) 

W ) 

(5.11) 


The  magnitude  and  phase  values  are  calculated  for  the 
frequency  response  for  the  transfer  functions  in  equation 
(5.8)  through  (5.11).  The  values  of  the  rigid  body  and  the 
flexible-fighter  aircraft  responses  are  practically  the  same. 
Table  8  illustrates  this  point. 


Table  8 


Frequency  Response 


G-COMMAND 

MANEUVER 

PITCH-POINTING  MANEUVER 

FREQ. 

PAD 

RIGID 

BODY 

FFA 

RIGID 

BODY 

FFA 

SEC" 

ldB 

ANGLE 

dB 

ANGLE 

dB 

ANGLE 

dB 

ANGLE 

0.2 

21.3 

-90.9 

21.3 

-90.8 

14.2 

-91.3 

14.2 

-91.2 

0.4 

15.3 

-92.0 

15.3 

-92.0 

8.23 

-92.8 

8.22 

-92.7 

0.8 

9.34 

-94.4 

9.34 

-94.4 

2.2 

-95.7 

2.2 

-95.7 

1.6 

3.42 

-99.8 

3.42 

-99.8 

-3.88 

-101.7 

-3.88 

-101.7 

2.2 

0.69 

-104.2 

0.69 

-104.2 

-6.76 

-106.2 

-6.76 

-106.2 

3.5 

-3.45 

-114.4 

-3.45 

-114.4 

-11.2 

-115.5 

-11.2 

-115.5 

5.0 

-7.0 

-126.1 

-7.0 

-126.1 

-14.9 

-125.1 

-14.9 

-125.1 

The  next  step  is  to  select  values  for  the  pilot  model  that 
allow  the  standard  of  performance  to  be  met. 

The  phase  angle  from  the  exponential  term  in  the  pilot 
model.  (57.3) (-0.3) («d)  ,  is  added  to  the  phase  values  of  0/F_ 
for  both  maneuvers.  This  establishes  the  open-loop  0/0 
magnitude  and  phase  information.  Figures  41  and  42  represent 
the  open-loop  amplitude-phase  plot  of  0/0_  for  the  g-command 
and  pitch-pointing  maneuvers,  respectively,  of  the  rigid  body 
and  the  flexible-fighter  aircraft  (FFA)  with  xp^  and  xp2  equal 
to  zero.  The  standard  of  performance  (closed-loop)  boundaries 
are  established  in  both  figures  on  the  Nichols  chart. 

First  examining  the  g-command  open-loop  plot  of  0/00 
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(Figure  41),  the  value  of  Kp  alone  can  not  meet  the  standard 

for  the  closed-loop  tracking  system  6 / 8 _ .  The  3.5  radians 

per  second  point  on  the  amplitude-phase  plot  requires  an 

increase  in  phase  angle,  indicating  lead  compensation  that 

must  be  added  by  the  pilot.  Figure  43  from  Reference  7 

provides  a  plot  of  amplitude-phase  curves  for  optimum  pilot 

compensation  which  simplifies  the  selection  of  xp^  and  xp2. 

Selecting  x  *0.8  and  x  =0.0  meets  both  the  bandwidth 
PI  Pi 

and  droop  requirements.  These  values  indicate  that  the 
pilot  must  generate  about  39  degrees  of  lead  at  a  frequency 
of  3.5  rad/sec.  The  resulting  (e/0  |  max  is  plotted  (as  a 
dot)  verses  the  pilot  compensations  angle  on  the  criterion 
graph  in  Figure  44.  The  point  falls  within  the  desired  PR 
boundary . 

Examining  the  pitch-pointing  open-loop  plot  of  e/e 
(Figure  42),  the  value  of  Kp  alone  can  not  meet  the  standard 
for  the  closed-loop  tracking  system  0/0_.  Again  lead 
compensation  by  the  pilot  is  needed  according  to  the 
amplitude-phase  plot.  Selecting  x  =  0.9  and  xp2  *  0.0 
meets  the  bandwidth  and  droop  requirement.  These  values 
indicate  that  the  pilot  must  generate  about  42  degrees  of 
lead  at  3.5  rad/sec  frequency.  This  point  is  also  plotted 
(as  a  triangle)  in  Figure  44.  According  to  the  point  on 
the  figure,  the  design  of  pitch-pointing  maneuver  would  not 
be  acceptable  to  the  pilot. 


Figure  43.  Amplitude-Phase  Curves  for  "Optimum"  Pilot 
Compensation 


Criterion  Graph  with  Design  Evaluation  Points 


The  closed-loop  resonance  gain,  | 0/ ©c |  yAX,  that  is 
plotted  in  Figure  44  comes  from  the  magnitude  of  Kp  and 
the  gain  established  by  the  (Tpls-t-l)  (tp2s+2) .  That  is 


K  *  K. 


r<TPl^2  *  1 

(Tp2U)2  +  1 


1  H 


(5.12) 


where  K  is  the  closed-loop  resonance  gain.  Therefore,  the 
gain  that  is  plotted  versus  the  pilot  compensation  angle 
comes  from  the  pilot  model  only.  If  the  gain  from  the 
pilot  model  can  be  reduced,  the  pitch-pointing  design  is 
acceptable  by  the  pilot.  A  suggestion  was  made  by  Mr. 

Finley  Barfield  of  the  Flight  Dynamics  Laboratory  at  Wright- 
Patterson  AFB,  Ohio,  to  modify  the  Neal-Smith  model  in 
Figure  38.  The  modified  model  is  shown  in  Figure  45. 
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„  _-0.3s  Tpls+1 

F  '  F 

s  —  s 

e 

e 

Y 

p  td2s+1 

m  g  •  • 

gain 

Fs 

_ 

Figure  45.  Modified  Neal-Smith  Model  for  Pitch 
Attitude  Control. 

This  modified  model  decreases  the  effort  required  by 
the  pilot  to  perform  a  pitch  rate  maneuver.  This  reduces 
the  tendency  to  produce  a  PIO  when  performing  a  longitudinal 
maneuver . 

Using  the  values  of  xpl  and  xp2  selected  previously 
for  the  pitch-pointing  open-loop  plot  of  0/0  ,  and  setting 
g  *  lO.OdB,  the  point  plotted  (as  a  square)  in  Figure  44 
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falls  within  the  desired  PR  boundary.  This  is  true  only  if 
the  gain  plotted  in  Figure  44  comes  from  the  pilot  model 
alone . 

If  the  new  model  is  a  fair  representation  of  the 
original  Neal-Smith  model,  the  only  Standard  that  is  needed 
to  be  met  is  that  the  pilot  compensation  angle  must  be 
between  -20  and  50  degrees.  The  reason  is  that  the  gain 
g  may  be  increased  or  decreased  to  ensure  that  the  gain 
required  from  the  pilot  model  falls  within  the  acceptable 
PR  boundary. 


VI .  Conclusion 


The  major  objective  of  this  thesis  has  been  to  deter¬ 
mine  how  a  flexible-fighter  aircraft  can  react  to  a  design 
that  is  originally  developed  for  a  rigid  body  aircraft. 

The  result  found  in  this  thesis  is  that  this  flexible- 
fighter  aircraft,  as  described  in  Chapter  II,  responds  the 
same  as  the  rigid  body  aircraft  when  using  the  same 
controller  design.  Even  when  the  flexible-fighter  model 
is  changed  through  a  range  of  dynamic  properties,  as 
described  in  Chapter  IV,  the  time  responses  are  nearly  the 
same  as  for  the  rigid  body  model. 

The  performance  obtained  from  the  digital  controller 
developed  with  Professor  Porter's  technique  is  very  good. 
Both  the  g-command  and  pitch-pointing  maneuvers  shown 
excellent  performance  with  their  individual  controllers. 
Even  when  a  time  delay  is  incorporated  to  represent  the 
controller  implementation,  both  the  rigid  body  and  the 
flexible-fighter  aircraft  simulation  time  responses  are 
not  significantly  affected  and  the  aircraft  performance 
is  very  satisfactory. 

According  to  the  Neal-Smith  Criterion,  the  g-command 
controller  is  a  pilot  acceptable  design.  Evaluating  the 
pitch-pointing  controller,  the  results  show  initially  that 
the  design  is  not  acceptable.  After  modifying  the  model 
(Figure  45),  the  pitch-pointing  controller  is  again 
evaluated  using  the  new  model  and  the  results  show  that  it 


is  a  pilot  acceptable  design. 

Equations  (5.8)  and  (5.9)  for  the  g-command  transfer 

function  0/F_  of  the  rigid  body  and  the  flexible-fighter 
s 

aircraft,  respectively,  have  an  unstable  real  pole  or 
complex  pole  pair.  Since  the  maneuver  only  lasts  8 
seconds,  the  unstable  poles  for  both  aircraft  do  not 
affect  the  time  response.  Because  the  poles  are  very 
close  to  the  origin,  the  pilot  is  able  to  control  the 
unstable  response.  Also,  a  g-command  maneuver  typically 
may  probably  not  last  more  than  8  seconds. 

The  transfer  function  6 /F  of  the  flexible-fighter 
aircraft  performing  the  pitch-pointing  maneuver  has  one 
unstable  pole  as  indicated  in  Equation  (5.11).  The 
maneuver  lasts  4  seconds  and  the  unstable  pole  does  not 
affect  the  response.  If  the  maneuver  lasts  longer  than 
10  seconds,  the  unstable  pole  can  probably  start  to  influence 
the  aircraft's  response,  but  the  pilot  can  still  be  able  to 
control  the  aircraft,  that  is,  to  break  off  the  maneuver. 
Usually,  as  in  a  dog  fight,  a  pitch-pointing  maneuver  may 
not  last  longer  than  10  seconds,  especially  if  the  aircraft 
is  in  a  turn  pulling  over  4  g's. 
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Appendix  A 
Subroutine  Fix 


Introduction 

The  purpose  of  the  Subroutine  Fix  is  to  make  the 
program  MULTI  user  friendly  when  inputing  matrices.  The 
way  this  works  is  that  the  user  will  no  longer  have  to 
input  the  matrix  from  beginning  to  end  if  an  error  is 
made.  All  the  user  has  to  do  now  is  to  give  the  row  and 
column  number  along  with  the  new  value.  This  is  explained 
in  more  detail  in  the  User's  Guide. 

User's  Guide 

The  options  which  ask  for  the  input  of  certain 
matrices  have  the  same  format  and  the  user  enters  the 
data  by  row  after  the  prompt  (>).  After  each  matrix  is 
entered,  its  values  are  automatically  echoed  back  for 
checking  (Ref.  9).  The  checking  tool  used  is  the  following 

IS  THIS  CORRECT _ YES,  NO,  $...> 

In  this  case  of  a  "$"  reply  the  option  is  terminated, 
although  the  option  is  terminated,  the  values  just  entered 
are  still  in  memory.  The  reply  "N  "  gives  the  following 
message  (the  user  can  also  enter  a  simple  "Y"  or  "N"): 

ENTER  THE  NUMBER  OF  CORRECTIONS  BEING  MADE...> 

The  number  of  corrections  correspond  to  the  number  of 
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wrong  values  in  the  matrix.  After  entering  the  number  of 
corrections,  another  message  appears  along  with  a  prompt 
(>). 


ENTER  THE  ROW,  COLUMN,  AND  NEW  VALUE  FOR  EACH  CORRECTION 
EX->  1, 2,-1. 5. . .PRESS  RETURN. 

> 

Only  one  correction  can  be  entered  at  a  time.  Once 
the  row,  column,  and  new  value  for  that  single  correction 
is  entered  the  user  must  press  the  return  key.  The  prompt 
( >)  appears  according  to  the  number  of  corrections  being 
made.  After  the  corrections  have  been  made,  the  matrix  is 
echoed  back  with  the  new  values  for  checking.  The  same 
checking  tool  is  used  again. 

Programmer 1 s  Guide 

The  Subroutine  Fix  is  used  to  make  MULTI  user  friendly. 
The  subroutine  is  located  in  the  Main  Executive  Overlay. 

The  following  fortran  statements  are  found  in  the  Subroutine 
Fix: 


SUBROUTINE  FIX.  (TS,IS,ID) 

DIMENSION  TS(10, 10) 

This  Initializes  the  subroutine  that  has  three  para¬ 
meters  that  are  described  below: 

TS  -  A  real  matrix  of  maximum  dimension  10  x  10 

TS  -  An  integer  value  denoting  the  #  of  rows  of  the 
matrix . 
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ID  -  An  integer  valued  denoting  the  #  of  columns  of  the 
matrix. 

These  parameters  are  transferred  from  the  previously 
entered  matrix. 


PRINT  '(/A)',  'ENTER  THE  NUMBER  OF  CORRECTIONS 
BEING  MADE . . . > ' 

READ*,  K 


First  the  statement  is  printed  out,  then  the  reply 
is  read  that  will  be  used  in  a  DO  Loop 


PRINT*,  'ENTER  THE  ROW,  COLUMN,  AND  NEW  VALUE  FOR 
EACH  CORRECTION. ' 

PRINT*,  'EX->1,2,-1.5  ...  PRESS  RETURN. ' 

DO  60  L=1 , K 
PRINT  '(/A)’,  '  >* 

60  READ* , I , J , TS ( I , J ) 

PRINT*,  '  ' 

END 


This  is  the  main  part  of  the  subroutine  which  actually 
corrects  the  errors  in  the  matrix.  The  Subroutine  Fix  has 
been  added  to  MULTI  throughout  the  program.  Everyplace  a 
matrix  is  inputted  the  Subroutine  Fix  has  been  added.  Below 
is  the  program  before  the  fix  was  added:  (this  is  not  an 
actual  matrix) 

1  PRINT*,  'ENTER  "Y"  MATRIX...  ,N, '  ROWS  WITH  ',K, 

'  ELEMENTS  EACH’ 

DO  2  I-1,N 

PRINT*,  'ROW  ',1,'  >’ 

2  READ*,  (Y(I,J),J-1,K) 
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+  CALL  ANSWER  (*1,*8001) 

CALL  MATPR ( Y , N , K ) 

PRINT* ,  'Y  MATRIX...' 

+  NOTE-  1st*  -  line  number  denoting  inccorect  data 
(Ref.  9)  directing  return  to  number  given.  (*1) 

2nd*  -  line  number  denoting  option  abort 
directing  return  to  end  of  option 
(*8001) 

This  is  how  the  program  looks  with  the  subroutine  daded. 

PRINT*,  'ENTER  "Y"  MATRIX. .. ;N, '  ROWS  WITH,',K, 

'  ELEMENTS  EACH' 

DO  2  II ,N 

PRINT*,  'ROW  ',1,'  >' 

2  READ*,  (Y( I , J) , J*1 ,K) 

GO  TO  5 

1  CALL  FIX  (Y,N,K) 

5  PRINT*,  'Y  MATRIX...' 

CALL  MATPR  (Y,N,K) 

CALL  ANSWER  (*1,*8001) 


Appendix  B 

Combination  of  States 


Introduction 

The  output  equation,  y«Cx.  can  be  considered  to  be 
the  combination  of  states  formed  by  the  "C"  matrix.  The 
output  for  a  MIMO  system  can  be  evaluated  from  the  state 
and  output  equations  (Ref  1.): 

• 

x  *  Ax  +  Bu  (B-l) 

£  *  Cx  (B-2 ) 

Taking  the  Laplace  transform  of  Equations  (B-l)  and  (B-2) 
yields 

sx(s)  -  x(0)  »  Ax(s)  +  BU(S)  (B-3) 

Y(s )  -  Cx(s)  (B-4) 

Solving  for  x(s)  in  Equation  (B-3) 

X(s)  -  £(s)x(0)  +  £(s)BU(s)  (B-5 ) 

where  $(s)  »  [sI^-A]  ”1.  Substituting  in  Equation  (B-5) 
into  Equation  (B-4)  yields 

Y(s)  -  C4(s)x(0)  +  C*(s)BU(s)  (B-6) 

The  output  equation  above  is  actually  the  combination  of 
states,  but  that  is  only  true  because  of  the  output 
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equation  Y(s)  =  Cx(s).  If  the  output  equation  happened 
to  be  Y(s)  =  Cx(s)  +  DU(s)  it  can  no  longer  be  just  the 
combination  of  states.  This  is  one  of  the  reasons  for 
developing  a  routine  that  takes  the  states  in  Equation 
(8-5)  and  forms  only  state  combinations.  This  is 
accomplished  by  forming  a  new  matrix  "Z"  and  multiplying 
it  by  the  states  determined  in  Equation  (B-5).  Another 
reason  for  this  routine  is  that  the  Program  MULTI  deals 
with  only  equal  number  of  inputs  and  outputs.  Therefore, 
until  now  only  the  combination  of  states  £  »  Cx  appearing 
in  the  output  could  be  observed. 

User's  Guide 

The  routine  for  combination  of  states  is  now 
contained  in  Options  31  and  32  for  terminal  plots,  and 
Options  34  and  35  for  calcomp  plots.  The  first  part  of 
the  User's  Guide  deals  with  Options  31  and  34.  The  second 
part  covers  the  use  of  Options  32  and  35. 

After  31  or  34  are  entered,  either  option  gives  the 
reply  (after  stating  whether  the  option  produces  a  terminal 
plot  or  a  calcomp  plot): 

PLEASE  CHOOSE  ONE  OF  THE  FOLLOWING: 

FOR  A  SINGLE  SAMPLING  TIME 

1.. .A  PLOT  OF  UP  TO  2  INPUT  AND  OUTPUT  PAIRS 

2.. . A  PLOT  OF  UP  TO  4  INPUTS  OR  OUTPUTS  OR  STATES 

3.. . A  PLOT  OF  UP  TO  4  DIFFERENT  SIMULATIONS 

(FOR  ANY  SINGLE  INPUT  OR  OUTPUT) 

OR  4... A  PLOT  OF  UP  TO  4  COMBINATION  OF  STATES 


ENTER  CHOICE  DESIRED  >4 


Choice  number  4  is  the  routine  that  enables  the  user 
to  form  conbinations  of  states.  When  choice  4  is  entered 
the  statement  below  is  obtained. 

CHOC IE  #4... YOU'VE  CHOSEN  TO  PLOT  COMBINATION  OF 
STATES 

After  the  statement  is  given,  the  following  question 
appears : 

ENTER  THE  NUMBER  OF  COMBINATIONS  OF  STATES > 

This  is  asking  for  the  number  of  different  combinations 
the  user  wants  to  plot.  Choice  number  4  permits  the  user  to 
enter  as  many  as  four  (no  more  than  four)  different  combina¬ 
tions  of  states.  When  the  number  of  combinations  is 
entered,  a  "Z"  matrix  must  be  formed  by  the  user.  This  "Z" 
matrix  must  contain  the  combinations  the  user  wants  to  see. 
This,  like  the  other  matrices,  has  the  same  format  and  the 
user  enters  the  data  by  row  following  the  prompt  (>).  Below 
is  an  example  of  a  5  state  model  in  which  the  user  wants  to 
add  states  2  and  4  together: 

ROW>0,l,0,i,0 

After  the  matrix  is  entered,  its  values  are  automatically 
echoed  back  for  checking.  The  chekcing  tool  used  is  the 
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same  one  described  in  Appendix  A. 

After  the  "Z"  matrix  is  entered,  the  following  reply 
is  given: 

FOR  GRID  ON  PLOT  ENTER  "0",  FOR  A  GRID  ENTER  "1"> 

This  reply  is  self  explanatory.  After  the  user 
decides  whether  or  not  a  grid  is  needed,  the  plot  for  the 
combinations  of  states  is  printed  out  on  the  terminal  for 
Option  31,  or  for  Option  34  a  calcomp  plot  is  developed. 

The  following  message  is  given  to  the  user  when 
Option  31  or  32  is  used: 

DO  YOU  WANT  A  LIST  OF  POINTS  USED  IN  PLOTTING? 

ENTER. . .YES  OR  NO. . .> 

In  case  of  a  "NO"  reply  the  option  is  terminated.  The 
reply  "YES"  gives  a  table  of  the  points  used  in  plotting. 

The  user  can  also  enter  a  simple  "Y"  or  "N".  The  following 
is  a  sample  format  of  a  table  produced  by  the  program. 

TIME  COMB  1 

0.000  0.  . 

.040  -.8704E-01 

.080  -.7281E-01 

.120  -.4229E-01 

The  maximum  number  of  points  used  in  plotting  is  one  hundred. 

Now  Options  32  and  35  deal  with  the  short  version  for 
plotting  either  on  the  terminal  or  developing  a  calcomp  plot. 
The  statement  below  is  given  for  either  option. 


COMB  2 
0. 

-.7894E-01 

-.6307E-01 

-.4558E-01 


ENTER  SHORT  VERSION  CHOICE. . .1,2, 3,4. .. > 

These  numbers  represent  the  same  numbers  found  in 
both  Options  31  and  34.  Choice  4  is  the  number  that  plots 
up  to  four  combinations  of  states.  From  there  on.  Options 
31,  32,  34,  and  35  have  the  same  format.  That  is,  the 
same  statements  are  given  as  before  according  to  the  option 
entered. 

Programmer's  Guide 

This  guide  does  not  cover  the  exact  location  of  the 
changes  made  or  the  additions  made  to  the  Program  MULTI . 

It  just  gives  the  changes  and  their  approximate  locations. 
The  numbers  that  were  changed  or  moved  are  not  mentioned. 
Therefore,  the  main  concern  of  this  guide  is  to  explain 
what  is  added  and  what  their  purpose  is.  The  reference 
numbers  given  (located  on  the  right  side  on  the  page)  are 
their  positions  relative  to  each  other  from  the  point  of 
location. 

In  the  overlay  that  contains  the  Options  30  thru  39 
the  following  changes  are  added. 

DIMENSION  Z(10,10)  3 

COMMON  /B  2  /  M  6 

These  statements  are  added  in  the  first  part  of  the 
program  OPTPLT.  The  statements  below  are  added  after  the 
comment  "SET  DATA  FOR  SHORT  VERSION  PLOT".  This  enables 
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The  user  to  use  the  short  version  to  plot. 


PRINT*,  'ENTER  SHORT  VERSION  CHOICE . . . 1 ,2 ,3 ,4 . . . > '  3 

IF  (ICODE.EQ.4)  GO  TO  1486  6 

From  the  comment  "SET  ICODE"  the  following  choice  is 
added  with  some  of  the  minor  changes  that  are  made. 

PRINT*,  'OR  4... A  PLOT  OF  UP  TO  4  COMBINATION  OF 
STATES ' 

IF  ( ICODE. LT.l. OR. ICODE. GT. 4)  GO  TO  1415 

GO  TO  (1425,1430,1435,1440)  ICODE 

Above  the  comment  statement  "SET  CHOICE,  SET  NUM, 
PRINT  MESSAGES  (1)"  the  following  are  added. 

1440  PRINT*,  'CHOICE  #4... YOU'VE  CHOSEN  TO  PLOT 
COMBINATION  OF  STATES* 

GO  TO  1486 

The  next  few  groups  of  fortran  statements  belong  to 
the  main  part  of  the  routine  that  combines  states.  These 
statements  are  located  three  lines  below  the  comment 
"OUTPUT(3)". 

1486  PRINT  '(/A)',  'ENTER  THE  NUMBER  OF  COMBINATIONS 

OF  STATES > ' 

READ*,  K 
PRINT* ,  '  ' 

The  first  two  statements  are  to  set  the  number  of 
combinations  of  states,  up  to  four,  that  can  be  performed. 
The  next  four  statements  are  used  to  generate  a  "Z"  matrix 


7 

10 

12 
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that  combines  the  states.  The  variable  "M"  is  a  common 


statement  in  which  the  number  of  states  in  the  system  is 
stored. 


PRINT*,  ‘ENTER  ”2'’  MATRIX. ..,  ‘K,  '  ROWS  WITH 
.  M, '  ELEMENTS  EACH' 

DO  1490  1-1  ,K 
.  PRINT*,  'ROW  ',1,'  >• 

1490  READ*,  (Z(I,J),J-1,M) 

PRINT* ,  ‘  ‘ 


The  next  statements  call  different  subroutines  and 
their  placement  in  the  routine  has  in  the  same  format 
shown  at  the  end  of  Appendix  A. 


GO  TO  190 

1495  CALL  FIX  (Z,K,M) 

190  PRINT*,  'COMBO  MATRIXZ. . . * 
CALL  MATPR  (Z,K,M) 

CALL  ANSWER  (*1945 , *8010) 


The  statements  below  are  the  actual  steps  used  in 
performing  the  combinations  of  states. 


DO  1510  I  -  1,N 
DO  1505  J  -  1,K 
S  -  0. 

DO  1500  L  -  1,M 
1500  S  -  S  +  X(I ,L+1)*Z(J,L) 
1505  PLMAT  ( I ,  J  )  -  S 
1510  CONTINUE 
ICLM  -  K 

CHOICE  -  '  COMBO' 


Figures  B-l  through  B-4  show  the  actual  set  ups  for 
each  matrix. 
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Figure  B-3.  Representation  of  the  multiplication  of  X(I,L+1) 
and  Z(J, L). 


Figure  B-A.  Representation  of  the  matrix  PLMAT(I,J). 
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Figure  B-4  is  the  format  of  matrix  PLMAT  in  order 
for  the  terminal  or  calcomp  plot  to  be  developed.  The 
statement  CHOICE  =  '  COMBO'  is  added  so  that  identifica¬ 
tion  of  the  combination  can  be  made  after  the  terminal 
plot  is  printed  out. 

The  next  statements  are  in  the  overlay  for  terminal 
plots.  One  common  statement  is  needed  along  with  two 
character  string.  These  two  statements  are  needed  if 
the  table  described  in  the  User's  Guide  is  going  to  be 
printed. 

COMMON  N/1B  13 C/  13C/  V(101,ll) ,X(101,11) 

CHARACTER  CHECK,  H(5)*6 

The  statements  below  are  needed  to  identify  the 
different  combinations  the  terminal  plot  illustrates. 

These  statements  are  located  after  the  comment  "EXIT  TO 
PRINT  STATEMENTS  WHEN  PLOT  IS  DISPLAYED". 

.  IF  (IC0DE.EQ.4)  GO  TO  1630 
1630  PRINT*,'  CURVE  ' , STRING (i : I ) , '  ABOVE  IS  ', 
choice, '  ' ,  I 

The  following  statements  are  used  to  print  the  table 
that  contains  the  data  for  plotting. 

1636  PRINT*  '  '  13 

PRINT*'  'DO  YOU  WANT  A  LIST  OF  POINTS  USED 
IN  PLOTTING?' 

PRINT  '(/A)',  'ENTER. . .YES  OR  NO. . .>'  15 
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READ  '(A)',  CHECK  16 

IF  (CHECK. EQ. 'N' )  GO  TO  1650  17 

PRINT*,  '  '  18 

H(l)  =  'COMB  1'  19 

H(2)  =  'COMB  2'  20 

H(3)  *  'COMB  3'  21 

H(4)  *  .'COMB  4'  22 

PRINT  5005,  'TIME',  (H ( I ) , 1=1 , ICLM)  23 

PRINT*,  '  '  24 

DO  1638  Irl ,N  25 

1638  PRINT  5010,  X(I,1),  (PLMAT(I ,K) ,K=1, ICLM)  26 

GO  TO  1650  27 

5005  FORMAT  (2X,A4,4X,4(A6,7X) )  28 

5010  FORMAT  (IX, F6 .3 , 1X,4(1E10 .4 ,3X)  29 

.  IF  (CHOICE. EQ. 'COMB' )  GO  TO  1636  33 

1650  I ENTRY  =0  34 


The  next  two  statements  are  in  the  overlay  for  the 
calcomp  plot.  Their  location  is  after  the  statement  SYS  = 
'SIM  ' / /STRING (J : J)//  ’ 

ELSEIF  (IC0DE.EQ.4)  THEN 
SYS  =  STRING (J : J)// '  ' 


127 


Appendix  C 

Forming  the  Measurement  Matrix 
And  Changing  the  Kq  and  Matrix 
For  the  Irregular  Method 


Introduction 

The  program  MULTI  is  used  for  designing  Digital 
Flight  Control  Systems  (DFCS).  This  is  accomplished  by 
inputting  the  plant  matrices  and  all  the  design  parameters 
(epsilon,  sigma  matrix,  alpha,  measurement  matrix,  etc.) 
into  MULTI,  calculating  the  Kg  and  matrices,  and  then 
running  the  simulation  option.  Option  26  uses  the  measure¬ 
ment  matrix  M  and  the  calculated  Kq  and  ^  matrices  to  run 
the  simulation  for  the  DFCS.  It  would  be  helpful  in  the 
design  to  be  able  to  change  one  of  the  design  matrices 
without  changing  the  others.  Option  14  did  allow  the  user 
to  change  the  measurement  matrix  M  but  it  would  then 
recalculate  the  Kq  and  K^  matrices.  Therefore,  Option  16 
is  created  to  input  just  the  measurement  matrix,  and  Option 
114  is  changed  to  allow  the  user  to  change  either  the  Kq  or 
K1  matrix,  or  both.  Option  14  for  the  irregular  method  is 
still  used  in  the  initial  calucation  of  the  Kq  and  K^ 
matrices. 

Option  118  is  changed  to  116  to  allow  the  user  to  see 
the  measurement  matrix  without  changing  it. 

User's  Guide 

Option  16,  which  is  used  to  input  the  measurement  matrix 
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M,  uses  the  same  format  as  the  other  matrices  and  the  user 
enters  the  data  by  row  after  the  prompt  (>).  After  the 
matrix  is  entered,  its  values  are  automatically  echoed 
back  for  checking.  The  checking  tool  used  is  the  same 
one  that  is  described  in  Appendix  A. 

Option  114  initially  prints  out  the  and  K1  matrices. 
Then  the  following  statements  appear: 

DO  YOU  WANT  TO  CHANGE  KO,  OR  Kl,  OR  BOTH? 

ENTER  1  FOR  KO,  2  FOR  Kl,  3  FOR  BOTH,  OR  4  FOR  NONE.> 

In  case  of  a  "4"  reply  the  option  is  terminated.  The 
"1"  reply  gives  the  following  statements: 

THIS  IS  TO  CORRECT  KO 

ENTER  THE  NUMBER  OF  CORRECTIONS  BENING  MADE> 

The  number  of  corrections  corresponds  to  the  number  of 
values  the  user  wants  to  change  in  the  matrix.  After 
entering  the  number,  another  message  appears  along  with  a 
prompt  ( > ) . 

ENTER  THE  ROW,  COLUMN,  AND  NEW  VALUE  FOR  EACH 
CORRECTION 

EX->1 , 2 , -1 . 5 . . . PRESS  RETURN. 

> 

Only  one  correction  can  be  entered  at  a  time.  Once 
the  row,  column,  and  new  value  for  that  single  correction 
is  entered,  the  user  must  press  the  return  key.  The  prompt 
(>)  appears  according  to  the  number  of  corrections  being 
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made.  After  the  corrections  have  been  made,  the 
following  statement  appears : 


NEW  KO  MATRIX... 


When  the  statement  appears  the  corrected  matrix  is  echoed 
back  for  checking.  The  checking  tool  used  is  the  same  as 
described  in  Appendix  A. 

The  "2"  reply  is  the  same  format  as  "1",  and  the 
reply  "3"  contains  both  and  . 

Programmer's  Guide 

The  fortran  statements  used  for  inputting  the  measure¬ 
ment  matrix  M  are  deleted  from  the  original  MULTI  program. 
The  following  fortran  statements  are  found  in  Option  16  for 
the  inputting  of  the  measurement  matrix: 


2016  PRINT  '(A/)',  'THIS  OPTION  FORMS  THE  MEASUREMENT 
MATRIX' 

PRINT* ,  '  ' 

PRINT*,  'ENTER  THE  M  MATRIX. ..' ,P, '  ROWS  WITH  ’ ,N-M, 
'  COLUMNS' 

.  DO  550  I  -  1,P 
550  PRINT*,  'ROW  ' ,1,  '>' 

GO  TO  560 

555  CALL  FIX  (MM,P,N-M) 

560  PRINT  '(/A/)’,  'MEASUREMENT  MATRIX... ' 

CALL  MATPR  (MM,P,N-M).  . 

CALL  ANSWER  (*555 ,*8002) 

FLAG  (16)  -  1 
GO  TO  8002 


The  I FLAG  (16)  -  1  is  an  error  message  used  to  deter¬ 
mine  if  the  matrix  has  been  inputted  when  Option  116  is 


used.  The  overlay  for  the  error  statements  is  changed 
in  two  places  to  the  following: 


3016  PRINT*,  'MEASUREMENT  MATRIX  MISSING...  SEE 
OPTION  #16  OR  #18' 

3018  PRINT*,  'ERROR  #18  STATEMENT* 


Option  114  is  changed  to  allow  the  user  to  change 
either  the  So  or  matrix  or  both.  The  following  state¬ 
ments  are  found  in  Option  114. 


C  OPTION  #114 

2114  IF  (METHOD. EQ. 'X' )  GOT  TO  1824 

IF  (METHOD . EQ . ' U ' )  PLTYPE  -  ' UNKNOWN ' 

I F  ( METHOD . EQ . ’ R ’ )  PLTYPE  *  ' REGULAR ' 

I F  ( METHOD . EQ . ' I ' )  PLTYPE  =  ' I RREGULAR ' 

I F  ( METHOD . EQ . ’ B ' )  PLTYPE  =  ' BSTAR ' 

PRINT*,  'CONTROL  MATRICES  ARE  FOR  PLANTS  WHICH  ARE', 
PLTYPE 

1824  PRINT  '(/A/)',  ’  KO  MATRIX...’ 

CALL  MATPR  (K0,M,P) 

IF  (IFLAG(lO).NE.l)  GO  TO  1825 
PRINT  '(/A/)',  '  K1  MATRIX...' 

CALL  MATPR  (K1,M,P) 

GO  TO  1827 

1825  IF  (ALPHA. EQ.l. AND. METHOD. NE. 'B')  THEN 

PRINT  '(1A/)',  *K1  MATRIX  IS  IDENTICAL  TO  KO  MATRIX' 
GO  TO  1872 
ENDIF 

PRINT  '(/A/)',  'Kl. MATRIX. .. ' 

CALL  MATPR  (K1,M,P) 

1827  PRINT*, 'DC  YOU  WANT  TO  CHANGE  KO,  OR  Kl,  OR  BOTH?’ 
PRINT*,  'ENTER  1  FOR  KO,  2  FOR  Kl,  3  FOR  BOTH,  OR 
4  FOR  NONE . > ' 

READ*,  CHECK 

IF  (CHECK. EQ. 4)  GO  TO  8015 
IFLAG(IO)  -  1 

IF  (CHECK. EQ. 2)  GO  TO  1828 
PRINT*,  'THIS  IS  TO  CORRECT  KO' 

25  CALL  FIX  (K0,M,P) 

PRINT*,  ’  NEW. KO  MATRIX...' 

CALL  MATPR  (K0,M,P) 

CALL  ANSWER  ( *25, *3015) 
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IF  (CHECK. EQ.l)  GO  TO  8015 
1828  PRINT*  , 'THIS  IS  TO  CORRECT  K1 ’ 
50  CALL  FIX  (K1,M,P) 

PRINT*, '  NEW  K1  MATRIX. . . ' 

CALL  MATPR  (K1,M,)). 

CALL  ANSWER  ( *50, *8015) 

GO  TO  8015 


The  IFLAG(IO)  *  1  statement  is  set  when  either  or  both  of 
the  Kq  and  matrices  are  changed.  This  determines,  if 
Option  114  is  used  again,  (if  alpha  -  1)  if  both  of  the 
matrices  should  be  printed  out  or  if  the  statement  "K1 
MATRIX  is  IDENTICAL  TO  K0  MATRIX’  should  be  used.  There¬ 
fore,  if  one  of  them  is  changed  both  matrices  are  printed 
when  Option  114  is  used.  The  IFLAG(IO)  is  reset  to  zero 
after  otpion  14  recalculates  the  Kq  and  matrices. 
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20.  aircraft.  The  same  design  is  then  used  on  a  flexible-fighter 
aircraft.  The  effects  of  the  controller  for  the  rigid  body  design 
when  used  on  the  flexible-fighter  are  evaluated  in  this  thesis. 

The  results  showed  excellent  performance  by  the  flexible-fighter 
aircraft. 

The  flexible-fighter  model  is  based  on  the  rigid  body  dynamics 
of  the  AFTI/F-16  combined  with  the  first  bending  mode  of  a 
flexible  aircraft.  The  bending  mode  of  the  flexible  aircraft  was 
selected  to  correspond  to  the  expected  characteristics  in  a  typica 
fighter  aircraft. 

Incorporating  a  time  delay  in  the  controller  to  represent  the 
expected  microprocessor  computational  time  delay  did  not 
significantly  affect  the  simulation  time  response  of  either  the 
rigid  body  or  the  flexible-fighter  aircraft. 

To  determine  the  acceptance  of  the  design  by  a  pilot,  the 
Neal-Smith  Criterion  is  used.  The  results  of  applying  this 
criterion  showed  that  the  pilot  would  accept  the  design  for 
both  the  rigid  body  and  the  flexible-fighter  aircraft. 
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